dates E os d 
C hInaA IV r1 


第 40 卷 第 2 期 KM 学 进展 Vol.40, No.2 
2022 4F 6 A PROGRESS IN ASTRONOMY June., 2022 


doi: 10.3969/j.issn.1000-8349.2022.02.01 


H 537 eo = 4E x 8473 A SR 


REH, BOR. KRA AR HIP 


(1. 中 国 科 学 院 云南 天 文 台 ， 昆 明 650216; 2. 中 国 科学 院 大 学 ， 北 京 100049; — 3. 中 国 科 学 院 紫金 山 
天 文 台 ， 南 京 210023; 4. 山东 大 学 空间 科学 研究 院 ， 威 海 264209; 5. 中国 科学 院 天 文大 科学 中 心 ， 
北京 100012) 


摘要 : 日冕 物质 抛射 (CME) 是 从 太阳 抛 出 的 大 尺度 磁化 等 离子 体 团 ， 携 带 大 量 磁 通 量 和 等 离 

体 进 入 行星 际 空间 。 进 入 行星 际 空间 的 CME 也 被 称 为 行星 际 CME (简称 ICME)， 当 它 朝向 地 
球 运 动 并 传播 至 地 球 附 近 时 ， 会 与 磁 层 相互 作用 产生 地 磁 暴 和 其 他 空间 天 气 现象 。 现 有 观测 结 
果 提 供 的 二 维 数据 不 能 完整 描述 CME 的 真实 磁场 结构 和 等 离子 体 分 布 。 为 了 更 好 地 预测 ICME 
到 达 地 球 的 时 间 以 及 对 地 球 和 周边 环境 可 能 产生 的 影响 ， 需 要 了 解 CME 的 三 维 结构 和 速度 的 3 
个 分 量 。 介 绍 了 根据 现 有 的 成 像 观 测 结果 对 CME 进行 三 维 重 构 的 方法 ， 包 括 基于 日 冕 仪 数据 和 
日 球 成 像 数 据 的 两 类 重 构 方 法 以 及 与 CME 成 像 重 构 关联 度 较 高 的 CME 驱动 激 波 的 三 维 重 构 方 
法 。 每 种 方法 适合 不 同 特点 的 CME 事件 ， 也 有 各 自 的 局 限 性 和 需要 满足 的 约束 条 件 。 我 们 对 比 
了 几 种 不 同 重 构 方 法 获得 的 结果 ， 发 现 这 些 方法 对 CME 的 传播 速度 和 传播 角度 的 估算 结果 都 比 
较 接 近 ， 说 明 这 些 方 法 的 可 靠 性 高 。 最 后 对 CME 三 维 重 构 的 热点 问题 和 发 展 趋势 进行 了 讨论 。 
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H 40 i HUT (coronal mass ejection, CME) 是 日 冕 大 尺度 磁场 结构 的 平衡 遭 到 破坏 
的 产物 ， 一 部 分 磁化 等 离子 体 (1016. g) 携带 大 量 磁 通 量 (107? ~ 10 Mx) 在 短 时 间 内 被 
高 速 向 外 抛 出 ， 并 释放 出 大 量 的 能 量 (107 J). CME 的 速度 一 般 在 200 ~ 3000 km/s 之 
间 ， 平 均 速度 大 约 是 500 km/s". CME 有 云 状 、 扁 状 和 泡 环 状 等 多 种 形状 。 其 中 最 典型 的 
CME 拥有 三 分 量 结 构 ， 包 括 一 个 明亮 的 前 沿 、 一 个 暗 腔 和 一 个 亮 核 ~。 图 1 为 太阳 和 日 球 
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KX (SOHO)" ”上 搭载 的 大 角度 分 光 日 园 仪 
(LASCO) ”观测 到 的 具有 “三 分 量 ” 结 构 的 
CME。 这 种 结构 经 常 被 观测 到 ， 但 是 目前 得 到 
的 可 见 光 波段 的 观测 结果 都 是 三 维 结构 在 二 维 
平面 的 投影 ， 而 从 不 同 角度 获得 的 投影 结果 是 不 
同 的 ， 所 以 目前 对 于 CME 三 分 量 的 真实 三 维 结 
构 还 不 是 很 清楚 。 由 于 CME 是 来 自 太 阳 的 大 尺 
度 磁 化 等 离子 体 结构 ， 在 行星 际 空间 中 传播 时 会 
对 周围 环境 产生 剧烈 的 扰动 。 当 这 些 庄 挟 着 等 离 
子 体 的 磁场 结构 到 达 地 球 时 ， 会 与 地 磁 层 相互 作 
用 ， 产 生 极光 ， 引 发 地 磁 暴 (目前 认为 最 严重 的 
2007/02/27 05:18 地 磁 扰 动 是 由 CME 引起 )， 对 地 球 周 围 环境 ( 即 
空间 天 气 ) 产生 扰动 ， 以 及 对 我 们 的 通信 系统 、 
电力 系统 、 高 空 飞机 ， 甚 至 输油管 道 造成 损害 ， 
同时 会 损坏 航天 设备 并 威胁 宇航 员 的 生命 安全 。 
因此 它 在 太阳 物理 以 及 空间 科学 领域 中 都 是 非常 重要 的 研究 课题 。 

由 于 地 面 观测 受 时 间 、 天 气 和 观测 波段 的 限制 ， 通 常 使 用 空间 设备 对 CME 进行 观测 。 
SOHO 卫星 ”是 迄今 为 止 观 测 到 CME 最 多 的 空间 观测 设备 ， 它 位 于 第 一 拉 格 朗 日 点 (L1) 
附近 ， 可 以 对 太阳 进行 24 h 不 间断 的 观测 。SOHO 所 携带 的 LASCO/C2 和 C3 HT aT 
以 对 太阳 周边 2.0 Ro ~ 30 Ro 范围 内 的 区 域 进行 观测 (这 里 Ro 代表 太阳 半径 )， 但 这 些 仪 
器 只 能 从 一 个 角度 对 CME 进行 观测 。2006 年 ， 日 地 关系 天 文 台 (STEREO)" 飞船 发 射 成 
功 ， 它 包括 A, B 两 艘 完全 一 样 的 飞船 和 载荷 ， 一 个 在 地 球 前 面 ， 一 个 在 地 球 后 面 绕 太阳 运 
动 ， 与 太阳 的 距离 分 别 是 0.9 AU 和 1.1 AU (这 个 距离 是 可 变 的 )， 其 实质 是 两 颗 人 造 行星 。 
STEREO-A 和 STEREO-B 的 分 离 角度 每 年 大 约 以 45° 增加 ， 其 中 所 携带 的 SECCHI™ 仪器 
包括 两 个 白光 日 园 仪 (COR1, COR2)、 一 个 极 紫 外 成 像 仪 (EUVI) 和 两 个 日 球 成 像 仪 (HI); 
主要 任务 是 从 不 同 角度 获得 CME 的 图 像 。 其 中 日 球 成 像 仪 可 以 对 内 日 球 层 的 ICME 进行 
成 像 观测 ”。 但 这 些 日 冕 仪 或 者 日 球 成 像 仪 都 只 能 获得 CME 三 维 结构 在 天 空 平面 的 二 维 投 
影 ， 不 能 完全 反映 CME 三 维 结构 ， 要 获得 CME 的 完整 几何 信息 ， 我 们 需要 根据 这 些 图 像 
和 相关 信息 ， 加 上 合理 假设 ， 还 原 出 CME 三 维 结构 。 这 一 过 程 被 称 为 CME 的 三 维 重 构 。 

根据 使 用 的 观测 数据 大 致 将 CME 成 像 重 构 方 法 分 为 两 类 : 基于 日 园 仪 观测 的 重 构 和 基 
于 日 球 成 像 仪 观测 的 重 构 。 这 些 方法 通过 使 用 不 同 的 观测 数据 进行 三 维 重 构 ， 得 到 的 三 维 参 
数 也 不 尽 相 同 。 基 于 日 园 仪 和 日 球 成 像 仪 的 观测 可 以 得 到 CME 的 三 维 形态 和 等 离子 体 分 布 
情况 。 我 们 将 在 第 2 章 详 细 介绍 不 同 的 CME 成 像 重 构 方 法 ， 在 第 3 章 简 单 介绍 两 种 CME 
驱动 激 波 的 重 构 方法 ， 在 第 4 章 对 其 中 一 些 CME 成 像 重 构 方 法 的 特点 和 适用 范围 进行 对 比 
和 深化 讨论 ， 然 后 在 第 5 章 对 本 工作 进行 总 结 ， 并 对 将 来 的 工作 进行 展望 。 


Hk: 图 片 来 自 SOHO/LASCO/C3。 
Eli CME 典型 “三 分 量 ” 结 构 
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2 CME 成 像 重 构 方法 


CME 具有 非常 复杂 的 动态 结构 ， 并 处 于 高 度 动态 的 运动 过 程 当中 ， 其 形态 和 大 小 各 有 


不 同 ， 有 膨胀 的 泡 环 状 、 
数 CME/ICME 都 表现 为 


“三 分 量 ”结构 、 扇 形 结构 等 。 在 行星 际 空间 的 传播 过 程 中 ， 大 多 
巨大 的 磁 通 量 绳 。 要 理解 并 最 终 预测 太阳 爆发 ， 就 必须 清楚 地 了 解 


CME 爆发 前 的 磁 通 量 强 内 部 结构 。 这 里 ， 人 磁 通 量 强 是 指 磁力 线 螺旋 缠绕 形成 的 磁场 结构 。 


样 的 结构 既 可 以 在 爆发 之 前 存在 ， 也 可 能 在 爆发 过 程 中 形成 。 一 般 认 为 爆发 前 存在 的 磁 


这 
通 量 绳 是 在 对 流 区 形成 ， 


然后 因 Parker 不 稳定 性 从 对 流 区 上 浮 至 光 球 和 色 球 ， 最 终 进入 日 


T 或 者 因 日 蝎 磁 结构 在 光 球 足 点 的 前 切 和 汇聚 运动 ， 驱 动 磁 重 联 在 足 点 发 生 ， 将 原先 简单 


- 


到 的 s 型 磁 结 构 等 ”"。 爆 


连结 的 磁 拱 转变 为 具有 螺旋 结构 的 复杂 磁场 位 形 ， 如 暗 条 (H) 和 爆发 前 经 常 可 以 观测 


发 过 程 中 形成 的 磁 通 量 强 ， 是 因为 剪 切 磁 拱 在 向 外 膨胀 的 过 程 中 ， 


其 后 方形 成 的 电流 片 中 发 生 的 磁 重 联 将 剪 切 磁 拱 转化 为 螺旋 缠绕 的 磁 结 构 所 致 。 在 实际 的 


ARM, B 


拱 和 磁 通 量 绳 完全 可 能 同时 存在 ”。 


无 论 哪 种 情况 ，CME 的 磁 结 构 本 质 上 都 是 一 个 大 尺度 的 磁 通 量 绳 ， 其 两 端 连 接 在 太阳 


上 ， 中 部 可 以 延伸 到 很 远 


的 行星 际 空间 ” (如 图 2 所 示 )。 


X 
电子 热流 密度 


PB 
地 球 


图 2 ”行星 际 日 园 物 质 抛射 在 行星 际 空间 传播 过 程 中 的 结构 示意 图 


CME/ICME 在 行星 


际 空间 中 传播 时 ， 随 着 与 太阳 距离 的 增加 而 膨胀 ;同时 


CME/ICME 在 传播 过 程 中 受到 太阳 风 或 者 其 他 结构 的 作用 还 会 发 生 偏转 和 形变 。Isavnin” 模 


拟 了 CME 磁 绳 的 完整 传播 过 程 及 其 在 太阳 风 中 发 生 的 变形 ， 表 明 CME 的 内 部 磁场 结构 在 
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其 传播 过 程 中 也 会 因 其 形变 而 发 生变 化 。CME/ICME 对 空间 环境 的 影响 与 其 在 特定 位 置 的 
磁场 结构 直接 相关 ， 而 磁场 结构 又 取决 于 CME 的 整体 三 维 几 何 和 形态 结构 。 因 此 ， 根 据 观 
测 数 据 对 CME 进行 正确 的 三 维 重 构 对 于 建立 可 靠 的 空间 天 气 预报 来 说 非常 重要 。 
2.1 BFA ROWAN Siz 
21.1 圆锥 模型 方法 

Zhao 等 人 "提出 了 一 个 圆锥 模型 (cone model) 来 估算 三 维 量 状 CME 的 几何 结构 和 运 
动 特性 。 在 圆锥 模型 中 ，CME 被 描述 为 流体 动力 学 等 离子 体 ， 其 特征 是 一 个 球形 的 自 相 似 
膨胀 结构 ， 即 具有 恒定 的 角 宽 、 传 播 方向 和 速度 。 尽 管 在 抛射 的 等 离子 体内 部 缺少 磁场 结 
构 ， 使 得 圆锥 形 模型 无 法 可 靠 地 预测 CME 在 地 球 或 其 他 空间 位 置 的 磁场 ， 但 是 这 个 模型 已 
经 成 功 地 用 于 预测 CME 是 否 到 达 地 球 以 及 到 达 地 球 的 时 刻 ”。 
圆锥 模型 基于 三 个 假设 : (1) CME 在 径 向 方向 上 以 几乎 恒定 的 角 宽度 通过 日 园 ; (2) 军 
JR CME 的 源 区 位 置 在 相应 的 活动 区 附近 ， 并 且 在 太阳 中 心 ，(3) CME 的 整体 速度 是 径 向 
的 ， 脱 上 胀 是 各 向 同性 的 。 另 外 圆锥 模型 选取 了 一 个 日 心 坐 标 系 (mn Yn Zn)» P 为 指向 地 
Bk, yn 指向 北 ， 而 on-yn 平面 定义 了 天 空 平 面 ， 再 建立 一 个 以 圆锥 顶点 为 中 心 的 右 圆锥 和 相 
应 的 坐标 系 (£o, Yc, zx)， 其 中 ze 沿 着 圆锥 的 对 称 轴 ， 而 yc-z 平面 平行 于 圆锥 的 底面 。 日 心 
坐标 系 中 的 (8, A) 定义 了 圆锥 的 方向 ，(%, A) 是 圆锥 中 轴 相 对 于 黄道 平面 的 经 度 和 纬度 角 ， 
而 圆锥 的 角 宽度 定义 为 2w (如 图 3 所 示 )。 


图 3 ”圆锥 模型 的 拓扑 结构 和 坐标 系 变换 


由 于 这 个 方法 需要 的 计算 资源 大 ， 且 模型 对 LASCO 图 像 的 拟 合 会 导致 结果 的 不 确 
定性 。 所 以 Xie 等 人 ”对 Zhao 等 人 的 圆锥 模型 进行 了 改进 ， 以 降低 计算 量 并 确定 晕 状 
CME 的 实际 速度 、 角 宽度 和 传播 方向 。 

为 了 更 好 地 从 观测 者 的 角度 (即日 心 坐标 系 中 ) 描述 圆锥 的 方向 和 传播 过 程 ， 需 要 将 
锥 坐标 系 变换 到 日 心 坐 标 系 。 为 此 ， 他 们 引入 了 一 个 过 小 坐标 系 (xt. yl, 21) (如 图 3 所 示 )。 
这 三 个 坐标 系 的 关系 如 下 : 首先 ， 将 轴 xy 和 加 逆 时 针 绕 轴 zn 转动 角度 a, VER zp 和 轴 oh 
EA, Wy 和 轴 伙 重合 ， 然 后 将 轴 2h (s) 和 轴 az 顺 时 针 绕 轴 ye (yl) 转动 角度 0 让 轴 z 与 
zc BA, fal 和 ze 重合 。 这 样 就 实现 了 日 心 坐标 系 到 圆锥 坐标 系 的 变换 ， 反 之 亦 然 ”。 


画 
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下 面 介绍 圆锥 模型 的 相关 参数 ， 在 锥 坐标 系 (me, y，xzc) 中 ，ze 与 圆锥 的 轴 重 合 ，7 是 贺 
锥 的 母线 ， 距 离 x 处 的 圆锥 截面 (如 图 3 中 的 圆 所 示 ) 可 表示 为 : 


Z,-—rcosu 
Ye=rsinwcosd , (1) 


Ze =rsinwsind 


其 中 ，6 = tan (z./y.) 为 圆锥 模 截 面 中 的 方位 角 。 将 式 (1) 转换 到 锥 模型 圆 截 面 在 平面 
x-y, ( 即 天 空 平 面 ) 上 的 投影 方程 : 


其 中 
h = r cos w cos 0 
a=rsinwsinĝð . (3) 
b=rsinw 
而 且 式 (3) 可 以 改写 为 : 
‘ _@ — enn a 
sinf = 7 >» 0 = sin (7) 
b 7 afta (4) 
tanw = (2) cos , w= tan (eos) 


R (2) 为 投影 在 天 空 平面 上 的 截面 的 椭圆 方程 ( 见 图 4), ZA a, b, h 分 别 为 其 短 轴 、 
长 轴 及 其 中 心 在 日 心 系统 中 的 原点 位 移 ， 可 通过 SOHO-LASCO-C2 和 C3 图 像 来 确定 。 由 
于 图 像 大 小 为 512 pixel 对 应 于 两 极 和 赤道 处 的 30 Ro, HA Xie 等 人 "将 C3 视 场 限制 在 
30 Ro 以 内 ， 而 且 a, b, h Alo 可 以 根据 图 4 所 示 的 几何 结构 计算 。 
传统 的 经 度 o 和 纬度 和 由 式 (5) 计算 : 


ee ( cos 0 sin a ) 


cos? 0 cos? a + sin? 0 


ó — tan! (SR 2) 


sin 0 


(5) 


[Ea 


此 ， 式 (4) 和 (5) 决定 了 圆锥 的 方向 与 角 宽 度 ， 图 5 给 出 了 利用 圆锥 模型 重 构 的 晕 状 
CME. 

如 果 au =b,h =0, Wi cos? = 0。 在 这 种 情况 下 ， 解 是 简 并 的 ， 因 为 角 宽 度 w 的 解 
不 是 唯一 的 。 当 圆锥 形 CME 的 方向 与 视线 平行 时 ， 就 会 发 生 这 种 情况 。 这 就 与 Zhao 等 
人 呈 遇 到 的 情况 相同 ， 这 时 Vaa 由 于 几何 原因 无 法 得 到 ， 所 以 Xie 等 人 ”使 用 经 验 公式 


C— 
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图 4 ”圆锥 模型 的 圆 形 截面 在 天 空 平面 上 的 椭圆 投影 


图 5 ”圆锥 模型 重 构 的 晕 状 CME (黑色 圆圈 ) 和 LASCO 在 2000 年 4 月 4 日 事件 中 的 图 像 的 对 比 


Via = 0.88Voxp ”来 弥补 这 一 不 足 。 
利用 式 (2) 和 (5) 得 到 CME 的 实际 径 向 速度 V, 与 投影 速度 Vi. 和 Vy, 之 间 的 关系 : 


EL re = C (cosa cos 0 — sinw sin 0 sin ô) 
dy. dr. (6) 
La = ap sw cos 


因此 


Chinay n "ZsEBBTII 
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dr Vy 
V= = i ] : 
dt | cosu cos — sin w sin Üsinó (7) 


同 理 ， 可 得 CME 实际 径 向 速度 V, 与 投影 速度 Von Vn 之 间 的 关系 为 : 


V. — dr Ven 

" dt coswcos 6 — sinw sina cos ô — sin@ cos a sin w sind (8) 
get U^ | 

p= 


di  cosw cos ð sin q&a + sin w cos q cos ô — sin f sin qa sin w sina 


这 里 Vons Vins Vax 和 Vy 是 天 空 平面 中 的 观测 速度 沿 轴 zh, yn, mL A y, 的 投影 速度 分 量 。 式 
(7) 和 (8) 都 可 以 计算 实际 径 向 速度 ， 通 过 比较 (7) 或 (8) 在 不 同 6 (PA) 下 的 结果 ， 可 以 检 
验 锥 模型 是 否 有 效 。 位 置 角 6 在 锥 坐标 系 下 定义 为 6= tan (ze/ye)， 极 角 (PA) 在 LASCO 
的 天 空 平 面 中 定义 为 P4 = —tan“'(a,/y,), PA 和 5 对 于 每 个 事件 都 有 一 个 复杂 而 独特 的 


— 


变换 关系 : 
PA = —tan™ (>) - 
Uh 
E cos 0 cos a cot w/ cos ô — sina — sin 0 cos a tan ô 
cos 0 sin a cot w/ cos 6 + cos a — sin 0 sin a tan ô 
不 同事 件 的 PA 与 方位 角 5 的 关系 取决 于 锥 的 三 个 自由 参数 ， 即 w, a AO. Xie ATH 
CME 的 实际 径 向 速度 ， 对 快速 尝 状 CME 和 慢 速 尝 状 CME 的 传播 时 间 分 别 进行 了 估算 ， 
结果 发 现 基于 模型 获得 的 结果 与 观测 结果 吻合 得 很 好 。 


圆锥 模型 方法 提高 了 建 模 的 精度 和 效率 。 这 是 首次 使 用 日 时 仪 数据 定量 地 确定 晕 状 
CME 的 实际 速度 、 角 宽度 和 源 区 位 置 。 这 些 参数 在 空间 天 气 模拟 中 至 关 重 要 。 锥 模型 在 使 
中 仍然 存在 一 些 不 确定 性 ， 比 如 CME 传播 时 间 的 预测 精度 依赖 于 CME 的 实际 速度 和 加 
速度 的 测量 ， 而 速度 由 观测 到 的 CME 在 不 同时 刻 的 位 置 推导 出 来 ， 加 速度 则 由 推导 出 来 
的 速度 进一步 推导 出 来 的 ， 另 外 ， 这 些 观测 与 推导 只 能 在 距离 日 心 30 Ro (LASCO/C3 视 
© 场 ) 的 范围 内 进行 测量 ，CME 出 了 这 个 距离 范围 之 后 就 无 法 测量 了 。 其 次 ， 圆 锥 模型 需要 的 
计算 资源 很 大 ， 且 模型 对 LASCO 图 像 的 拟 合 会 导致 结果 的 不 确定 性 ，Xie EA" 改进 了 
此 方法 ， 降 低 了 计算 量 并 确定 了 CME 的 实际 速度 、 角 速度 和 传播 方向 。 
2.1.2 ”偏振 比方 法 
偏振 比 (polarization ratio, PR) 方法 是 Moran 和 Davila ™ 首先 提出 ， 将 LASCO-C2 
的 偏振 观测 数据 转化 为 与 天 空 平面 (POS) 的 三 维 距离 。 后 来 Dere 5E A 利用 高 时 间 分 辩 
率 (1h) 的 LASCO 偏振 测量 对 PR 方法 进行 了 验证 ， 获 得 了 清晰 的 CME 结构 。Mierla 等 
AP" fii Moran 等 人 ”将 该 方法 应 用 于 STEREO 日 冕 仪 的 偏振 观测 ， 证 明 利 用 偏振 效应 可 
以 重 构 三 维 CME. Dai 4$ A" 利用 PR 方法 计算 了 CME 的 三 维 质量 ， 改 进 了 传统 CME 
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质量 的 计算 方法 ， 获 得 了 更 精确 的 CME 质量 。Lu 等 人 ”将 PR 方法 与 GCS 方法 相 结合 ， 
提出 了 一 种 新 的 GCS-PR 方法 来 获取 从 地 球 单一 角度 观测 到 的 CME 三 维 参数 。 

| 园 的 白光 辐射 来 自 自由 电子 对 光 球 辐射 的 散射 ， 因 此 ， 日 园 的 光 信号 是 偏振 的 。 我 们 
可 以 从 具有 不 同 偏振 态 的 辐射 强度 比 中 估算 出 有 效 散 射 角 ， 再 利用 该 散射 角 计 算出 与 天 空 
平面 的 距离 。PR 方法 的 物理 原理 是 汤姆 孙 散 射 ， 其 散射 截面 取决 于 散射 方向 与 电场 天 量 之 
间 的 夹 角 。 从 光 球 发 出 的 光 是 没有 偏振 的 ， 它 可 以 分 解 成 两 个 相等 且 相 互 垂 直 的 分 量 ， 一 个 
垂直 于 散射 面 ， 强 度 为 Kan， 另 一 个 平行 于 散射 面 ， 强 度 为 Raa。 根 据 汤 姆 孙 散 射 理论 ，TFian 
与 散射 角 x FEI, Iaa 正比 于 sin? xe MIRER wn、 总 亮度 ho 和 偏振 度 已 定义 为 ”: 


pol = Tian m Trad , (9) 
Liot = Lian + Trad = 2I = Ipa E (10) 
I ol 
pus ; 11 
Liot ( ) 
TO. = 
Tan = 5 nf dzNe(p, z)|(1— u)C +uD|, (12) 
0 
TO i ND 
ina = 5 Io dzNe(p, z)(1 — u)A + uB]sin^x , (13) 
0 


其 中 ， 是 太阳 圆 面 中 心 的 强度 ，o。 RAN LIRE, udi ESI AB (Dai 等 人 ” Hy 
u = 0.63), N 是 电子 密度 ，p 是 沿 天 空 平面 的 投影 距离 LA 6), 2 是 散射 点 沿 视 线 方向 
与 天 空 平面 的 距离 ，w 是 散射 点 在 太阳 圆 盘 上 所 张 的 半角 ，A, B, C, D 是 van de Hulst A 


x. 
A(r) = cos Q(r) sin? Q(r) ， (14) 
B()- -4 {1 3sin? 00) ER x pressit Gr] xm [EERON | as 
a 3 ae ose) (16) 
O Din) = -g {5+ sin? or) - EO- sin? oo] x m ELO), am 


HF, # OQ A sin Q(r) 2 1/r 给 出 (7 是 以 太阳 半径 Ro 为 单位 的 日 心 距离 ，7r? = p? + 27). 
对 于 每 一 条 视线 ，PR 方法 假设 所 有 贡献 于 和 an M Da 的 电子 都 位 于 同一 个 位 置 (po, z0) (如 
图 6 中 的 Pl 所 示 )， 称 之 为 等 效 散 射 中 心 ; 那么 在 等 效 散 射 中 心 的 电子 密度 为 N。(po, zo) 
而 且 偏 振 度 P 可 表示 为 偏振 亮度 与 总 亮度 的 比值 : 


[1 +u)A + uB] sin? x 
2[(1 +u)C + uD] — ((1— u)A +uB]sin? x ' 


(18) 
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Hh, A, B, C, DR x WE po, zo 的 函数 ，P 和 po 可 以 从 观测 中 得 到 ， 求 解 式 (18) 可 以 
得 到 散射 点 在 视线 方向 上 与 天 空 平面 的 距离 zx。 这 个 方法 的 缺点 是 得 到 的 结果 只 是 p. 55 
效 散 射 中 心 与 天 空 平面 距离 的 符号 具有 不 确定 性 ， 即 等 效 散 射 中 心 可 能 是 在 太阳 的 正面 也 
可 能 是 在 背面 。 我 们 可 以 利用 另外 一 个 视点 的 观测 结果 四 ， 或 者 利用 CME 源 区 的 EUV 观 
测 图 像 来 确定 zo 的 符号 。PR 重 构 不 能 得 到 关于 CME 在 视线 方向 上 的 分 布 ， 提 供 的 信息 只 
是 沿 视线 方向 上 的 CME 等 离子 体 密度 的 加 权 平 均 四 。 


ap 


Obs 


LE 


E: x 是 P1 处 的 散射 角 ，z 是 散射 点 沿 视线 方向 与 天 空 平面 的 距离 ，p 是 沿 天 空 平面 的 投影 距离 ，r 是 日 心 


Dr 


W 


6 ” 沿 视 线 方向 的 汤姆 孙 散 射 几何 示意 图 


~ 


Mierla 等 人 “利用 PR 方法 重 构 了 2007 年 8 月 31 日 的 CME， 如 图 7 所 示 : a) 和 b) 
表示 PR 方法 应 用 STEREO-A/COR1 数据 的 重 构 结果 ，c) Md) 则 表示 应 用 STEREO- 
B/COR1 数据 的 重 构 结果 ; a) 和 c) 是 侧 视 图 ，p) 和 d) 则 是 正视 图 。 从 图 7 可 以 看 出 重 构 
出 的 CME 只 有 一 个 粗糙 的 外 表面 ， 无 法 提供 任何 关于 CME 深度 层次 结构 的 信息 。 不 过 这 
种 方法 的 优点 是 只 需要 单一 角度 的 观测 资料 就 可 以 进行 重 构 。 

2.1.3 优化 圆柱 壳 方 法 

Thernisien 等 人 ”提出 的 一 种 磁 通 量 绳 正 演 模拟 技术 ， 称 为 优化 圆柱 过 方法 (graduated 
cylindrical shell, GCS)， 也 称 为 羊角 型 方法 ， 是 在 冰激凌 方法 的 基础 上 发 展 起 来 的 。 目 前 
冰激凌 方法 主要 应 用 在 曼 状 CME 的 三 维 重 构 ， 而 GCS 方法 假设 CME 的 膨胀 是 自 相 似 的 ， 

源 区 的 中 性 线 和 CME 的 角 宽度 基本 不 变 ， 对 CME 的 形态 进行 拟 合 反 演 。Cremades 和 
Bothmer (CB04)” 采用 锥 模型 对 124 个 CME 磁 通 量 绳 进行 了 系统 研究 ， 分 析 了 其 源 区 特 
征 与 LASCO 观测 到 的 CME 形态 之 间 的 关系 ， 并 推导 出 CME 的 投影 效果 。Thernisien 等 
AP" 基于 CB04 结果 对 CME 的 三 维 形态 进行 重 构 ， 获 得 了 OME 的 传播 方向 和 三 维 传播 
速度 。 

GCS 模型 中 的 CME 形态 类 似 于 羊角 管 ， 由 两 部 分 组 成 ， 即 磁 通 量 管 本 身 和 支撑 它 的 
两 个 锥 形 腿 ， 腿 的 两 端 连接 到 太阳 表面 。 图 8 是 GCS 模型 用 到 的 磁 通 量 强 结构 及 其 空间 位 
置 示意 图 ” ， 它 有 9 个 参数 ， 其 中 6 个 是 自由 参数 ， 分 别 是 : (1) CME 的 倾角 7， 即 相应 的 
瞳 条 或 者 磁 中 性 线 在 日 面 上 与 相 邻 纬度 线 之 间 的 夹 角 : (2) CME 的 半角 宽度 c， 与 暗 条 或 
者 磁 中 性 线 的 长 度 相 关 ， 当 a =O, GCS 模型 与 “冰激凌 ”模型 相同 ; (3) 磁 通 量 管 的 高 
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c) d) 


VE: a c) 为 重 构 结 果 的 侧 视 图 ，b), d) 为 正视 图 。a) 和 b) 表示 PR 方法 应 用 STEREO-A/COR1 数据 的 重 
构 结果 ，c) 和 d) 表示 应 用 STEREO-B/COR1 数据 的 重 构 结果 。 


图 7 ”PR 方法 的 重 构 结果 


度 即 腿 的 高 度 h; (4) 表征 磁 通 量 管 粗细 的 参数 x， 定义 为 磁 通 量 管 表面 的 一 点 到 轴线 的 垂 
直 距 离 与 其 到 日 心 的 距离 比 ， 即 « = sin 6; (5) Carrington 坐标 系 下 的 CME 源 区 的 经 度 o; 
(6) Carrington 坐标 系 下 的 CME 源 区 的 纬度 4。 这 些 参数 中 ， 前 四 个 定义 了 磁 通 量 绳 的 几 
何 结构 ， 最 后 两 个 确定 了 磁 通 量 绳 的 源 区 位 置 。 有 关 模 型 几何 结构 的 具体 说 明 见 文献 台 . 
通过 对 上 述 参 数 的 初始 化 ， 首 先 获 得 CME 在 天 空 平面 的 投影 ， 得 到 CME 的 三 维 结 
构 。 获 得 形态 拟 合 之 后 ， 再 选择 高 斯 拟 合 函 数 ” 对 CME 边缘 电子 密度 分 布 进行 假设 ， 即 : 


N,(d) - Nel CR], (19) 


N 
4 
F 


Otailing(d < a 
n=] tsiing 4 < a) (20) 


Cleadling(d 之 a) 


RF, d 是 壳 层 内 外 任意 点 到 骨架 (图 8 a) 中 的 虚线 ) 距离 。 该 函数 带 来 了 拟 合 需要 的 另 
外 3 个 参数 : 电子 密度 Ne > 高 斯 特征 长 度 Otailing All Oleadling。 该 函数 可 以 是 非 对 称 高 斯 分 
布 ， 电 子 密度 主要 位 于 训导 上 ， 这 种 不 对 称 性 会 使 模拟 结果 更 符合 实际 情况 。Thernisien 等 
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a) b) c) 
图 8 ”GCS 模型 中 磁 结构 的 正面 (图 a))、 侧 面 (图 b)) 和 空间 位 置 (图 c)) RED 


全 用 此 方法 重 构 了 26 个 CME， 并 与 Colaninno 和 Vourlidas” 确定 的 CME 方向 做 比较 ， 
除了 一 个 晕 状 CME 外 ， 其 他 结果 都 比较 相似 ， 传 播 方向 的 差别 在 10° 以 内 。 

图 9 展示 的 是 2008 年 4 月 26 日 的 一 个 事件 。 在 STEREO-A ( 右 ) 的 视 场 中 CME 向 东 
传播 ， 而 在 STEREO-B ( 左 ) 的 视 场 中 它 以 一 个 晕 状 CME 的 面目 出 现 。 由 于 STEREO-A 
和 STEREO-B 两 个 日 晃 仪 的 观测 角度 不 同 ，CME 前 沿 的 表现 也 不 同 ， 所 以 他 们 采用 了 手 
绘 前 沿 来 拟 合 GCS 模型 (图 9 第 二 排 )， 这 个 例子 的 拟 合 优 度 高 达 83 欠 。 然 而 对 于 另外 一 种 
情况 ， 如 图 10 展示 的 2007 年 12 月 31 日 爆发 的 三 个 CME， 尽 管 Thernisien 4$ AP* 已 经 将 
前 沿 分 为 两 个 部 分 来 进行 拟 合 ， 但 由 于 CME 的 前 沿 被 强烈 扭曲 ， 效 果 仍 然 很 差 ， 最 低 拟 合 
TERA 3496. 

GCS 模型 需要 两 个 不 同 视点 的 图 像 ， 最 适合 的 观测 设备 就 是 STEREO 卫星 ， 但 由 了 
2014 年 STEREO-B 卫星 不 工作 了 ， 所 以 目前 在 使 用 GCS 模型 时 一 般 同时 采用 STEREO-A 
和 SOHO/LASCO 的 数据 来 进行 拟 合 。 这 个 模型 的 最 大 优点 是 可 以 提供 CME 的 三 维 初 始 
速度 和 初始 方向 。 在 此 基础 上 Zhao 5 A" 采用 Liu 4$ A 提出 的 简单 解析 模型 ， 使 得 对 
CME 的 研究 不 再 局 限于 1 AU 以 内 ， 得 到 了 从 太阳 附近 到 5.34 AU 处 的 距离 和 时 间 的 关系 。 
当然 GCS 也 有 不 足 之 处 : 它 假设 CME 都 是 羊角 形 的 。 很 显然 ， 在 实际 情况 中 ， 很 多 CME 
都 不 是 羊角 形 的 。 对 2007 年 12 月 31 日 事件 的 拟 合 就 能 说 明 一 些 问题 。 另 外 ，GCS 只 能 针 
对 单个 CME 进行 模拟 ， 难 以 应 付 短 时 间 内 连续 爆发 的 多 个 CME。 对 连续 爆发 的 CME 的 
拟 合 结果 不 理想 ， 不 能 很 好 地 区 分 不 同 的 CME™。 

2.1.4 三 维 日 网 磁 绳 殷 射 方法 


~ 


Most] 等 人 ™ 对 太阳 风 


FP 磁 云 的 重 构 也 广泛 应 用 了 GS 方法 ， 


并 对 磁 绳 的 重建 提出 了 


另外 一 种 正 向 建 模 的 3 


Rope Ejection, 3DCORE 
除了 可 以 重 构 CME 还 可 以 得 到 地 磁 指数 (Dst)。 Fi 


和 经验 模型 ， 称 为 三 维 日 


MM AL Age 


ae XX XE UB 


EE) 


hA te AY (Three-Dimensional Coronal 


一 个 包含 行星 际 传播 以 及 CME 减速 和 膨胀 的 模型 ， 


T 


MARAT, CME 中 的 


向 磁场 在 


传播 过 程 中 大 致 守恒 ， 可 以 提供 足够 的 南 向 磁场 分 量 ， 与 地 球 磁 层 顶 中 的 北向 磁场 分 量 发 生 
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ik: 第 一 排 为 STEREO-B (£) 和 STEREO-A (4) 的 观测 图 像 ， 第 二 排 展 示 了 在 观测 结果 上 手绘 出 的 前 沿 ， 
第 三 排 为 利用 GCS 方法 重 构 获 得 的 结果 。 


9 “利用 GCS 技术 对 2008 年 4 月 26 日 爆发 的 CME 进行 重 构 的 结果 四 


Ae 
E 


jt: 第 一 排 为 STEREO-B (Jz) 和 STEREO-A (4i) 的 观测 图 像 ， 第 二 排 为 利用 GCS 方法 重 构 获 得 的 结果 。 


10 ”利用 GCS 技术 对 2007 年 12 月 31 日 爆发 的 CME 进行 重 构 的 结果 
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人 磁 重 联 之 后 ， 将 CME 带 来 的 磁 通 量 和 能 量 送 入 磁 尾 ， 引 发 强烈 的 地 磁 暴 。 所 以 ， 根 据 重 构 


结果 给 出 的 CME 轴 向 场 强度 ， 就 可 以 推导 得 出 CME 到 达 地 球 附 近 之 后 其 中 的 南 向 磁场 强 


度 ， 从 而 估算 出 地 磁 暴 的 强度 和 相应 的 Dst。 
3DCORE 模型 中 的 CME 几何 形状 如 图 11 
所 示 ， 它 有 2.5-D Gold-Hoyle"^ “模型 中 的 磁 通 
量 绳 结构 ， 横 截面 为 圆 形 ， 其 磁场 具有 均匀 扭曲 
的 位 形 ， 磁 通 量 强 两 端 连接 在 太阳 上 ， 截 面 和 整 
体形 状 都 是 圆 的 ， 但 横 截 面 的 半径 是 变化 的 ， 接 
到 太阳 两 端的 截面 半径 逐渐 趋 于 零 。 对 实际 事 
件 进 行 重 构 时 ，CME 的 初始 速度 和 传播 方向 由 
STEREO 的 HI (日 球 成 像 仪 ) 或 者 COR (白光 
HZA) 的 观测 结果 给 出 ， 磁 通 量 绳 的 手 征 性 和 
向 场 的 方向 由 CME 源 区 磁 图 和 极 紫 外 单 色 像 
得 出 。 其 前 沿 的 运动 遵循 drag-based Eug" ™ 
(简称 DBM 模型 ) 给 出 的 规律 ， 其 他 部 分 以 自 相 
似 的 方式 膨胀 ， 而 磁 通 量 强 中 轴 上 的 磁场 随 距 
离 的 增加 而 下 降 。 
3DCORE 模型 主要 用 于 研究 CME 在 行星 
际 间 的 传播 ， 以 及 预测 CME 在 1 AU 处 的 速 


: 


iE: 由 一 个 连接 到 太阳 的 锥 形 圆 环 (橙色 ) 

组 成 ， 黄 点 代表 太阳 ， 绿 点 代表 地 球 。 

11 3DCORE 模型 中 CME 磁场 结构 的 几何 
形状 


度 和 磁场 信息 。 它 的 初始 参数 ， 如 磁 通 量 绳 的 初始 位 置 和 CME 初始 速度 ， 可 以 根据 


STEREO/HI 观测 结果 获得 ， 或 者 由 GCS 建 模 结果 得 到 。 磁 场 的 缠绕 方式 可 直接 设 定 ”， 


或 者 由 观测 得 到 ”。 最 主要 的 一 个 参数 是 DBM 


的 阻力 参数 ， 其 范围 为 0.05 ~ 2"* 5, 3x 


个 参数 对 于 CME 到 达 时 间 的 预测 非常 重要 。 由 3DCORE 方法 得 到 的 结果 可 以 推断 出 磁 通 
量 绳 中 磁场 的 南 向 分 量 ， 这 对 预测 地 磁 暴 至 关 重要 ”。 但 是 3DCORE 方法 对 CME 前 沿 的 


运动 采取 DBM 模型 ， 而 阻力 参数 的 设 定 则 是 根 
CME 的 实际 情况 有 偏差 。 
2.1.5 三 角 测量 方法 


据 经 验 所 得 ， 往 往 导 致 模型 给 出 的 结果 与 


从 日 狗 仪 观测 数据 推算 CME 的 速度 会 受到 投影 效应 的 影响 ， 其 大 小 取决 于 CME 源 区 
在 太阳 表面 的 位 置 。 应 用 三 角 测量 的 方法 可 以 确定 CME 的 源 区 位 置 ， 消 除 投影 效应 ， 从 而 


获得 真实 速度 。Inhester 详细 描述 了 以 两 点 观测 确定 CME 的 三 维 空间 位 置 的 三 角 分 析 方 


沿 来 重 构 CME 的 空间 位 置 。 


ik. Temmer 等 人 "利用 STEREO 和 LASCO 两 个 位 置 的 观测 结果 ， 通 过 跟踪 CME 的 前 


三 角 测 量 法 的 基本 思路 如 图 12 所 示 。 相 关 参 数 包括 两 个 航天 器 各 自 看 到 的 天 空 平面 之 
间 的 夹 角 4, CME 前 沿 在 STEREO-A 和 LASCO 各 自 的 天 空 平面 上 的 投影 与 太阳 的 距离 
分 别 为 doa 和 dor» CME 前 沿 与 太阳 的 真实 距离 为 4，CME 在 LASCO 和 STEREO 获得 
的 图 像 中 的 延展 角 分 别 为 入 和 a。 这 里 延展 角 是 飞船 和 CME 前 沿 之 间 连 线 与 该 连 线 在 赤道 


面 上 的 投影 之 间 的 夹 角 ， 而 且 延 展 角 可 以 转换 为 


E 离 ， 比 如 对 于 LASCO 的 观测 dor = 216A 
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Ro ””。 根 据 这 个 关系 ， 我 们 可 以 先 获得 dor 的 大 小 ， 然 后 计算 STEREO-A 5 CME 前 沿 
的 延展 角 aw。 由 图 12 可 知 ，A。 是 CME 前 沿 上 的 一 点 在 STEREO-A 的 天 空 平面 上 的 投影 ， 
Al 是 A, 在 赤道 面 上 的 投影 ，CME 前 沿 上 的 同一 个 点 在 LASCO 天 空 平面 的 投影 为 L。.， 而 
L 则 是 Le 在 赤道 面 上 的 投影 ;点 万 则 是 直线 LL. 与 直线 AA, 的 交点 ，D 在 赤道 面 上 的 投 
WAC; 日 心 位 于 O 点 。 显 然 点 C 也 是 直线 LL 与 直线 AA 的 交点 。 角 AAA =a, ff 
LLL! = à; € J LAL 5 LL, ZIMA, « 为 OC 5 LC 之 间 的 夹 角 。 由 于 LOC 构成 一 
个 平面 三 角形 ， 因 此 +r +o = 180°% CME 源 区 在 日 面 上 的 经 度 和 纬度 分 别 用 vo 和 6 表 
示 ， 根 据 这 样 的 设置 和 安排 ， 我 们 可 以 得 到 如 下 结果 : 


d 
€ = arctan (=) ; (21) 
TAU 
7=180-e-y , (22) 
= 23 
加 i. sine , (23) 
> ô = arctan [tan(PA — A)sin(y+ya)] , (24) 
co 
AS p? — d? +riy —2sindcosdcos(p+ ya) , (25) 
2 2 d? 
(Arr) (26) 
2rau /p 


通过 改变 源 区 的 经 度 PAE 0. FEILER, O EAR BUE 5 a BIER RB ES 2L 
间 的 最 小 偏差 ， 从 而 得 到 对 CME 源 区 位 置 的 最 佳 估计 值 。 利 用 测量 的 投影 距离 (doA)， 将 
导出 的 经 度 插入 式 (21) — (23) 中 ， 就 可 以 估算 出 真实 距离 4。 需 要 注意 的 是 三 角 测 量 方法 
两 个 假设 ， 一 个 是 CME 的 速度 不 受 CME 膨胀 的 影响 ， 另 一 个 假设 是 CME 党 径 癌 传播 。 
很 显然 ， 这 些 假 设 在 有 些 情况 下 并 不 合理 ， 这 也 将 导致 相应 的 结果 出 现 较 大 的 偏差 。 另 外 由 
于 CME 图 像 是 沿 不 同 视 线 的 亮度 积分 的 结果 ， 所 以 在 两 幅 CME 图 像 中 的 对 应 点 的 识别 较 
困难 ， 与 实际 情况 有 偏差 ， 导 致 无 法 真正 实现 CME 的 三 维 重 构 ， 只 能 获取 CME 的 三 维 传 
秋 方 向。 不 过 使 用 三 角 测 量 方法 可 以 量化 和 纠正 大 多 数 测量 中 的 投影 效应 ， 而 不 需要 事先 得 
c 知 CME HZR. E Temmer 等 人 研究 的 11 个 事件 中 ， 大 多 数 事 件 的 重 构 结 果 都 
© 很 好 。 如 果 对 CME 的 距离 测量 没有 问题 ， 并 且 可 以 清楚 地 辨认 出 CME 在 三 个 仪器 中 相同 
的 结构 ， 那 么 三 角 测量 法 是 比较 有 效 的 。 目 前 三 角 测量 在 快速 佑 计 源 区 、 方 向 和 修正 投影 效 
应 中 都 取得 了 很 好 的 结果 。 
2.1.6 ŠAH k 
CME 重 构 中 还 有 一 种 常用 的 方法 是 点 共 斩 (tie point, TP) 方法 。TP 方法 需要 使 用 
STEREO 双星 提供 的 观测 数据 来 进行 三 维 重 建 。 两 个 STEREO 航天 器 和 日 田中 的 一 个 点 
定义 了 一 个 对 极 面 ” ， 当 目标 在 STEREO-A 和 STEREO-B 的 图 像 中 同时 出 现时 ， 就 可 以 
三 角 测 量 法 来 获得 其 三 维 日 心 坐标 系 中 的 位 置 ， 如 图 13 所 示 。 两 束 光 线 在 目标 的 位 置 相 
交 ， 寻 找 对 应 点 可 以 通过 眼睛 来 完成 ” (最 多 的 是 跟踪 前 沿 的 显著 特征 ， 如 密度 变化 、 形 状 


T 


cm 
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pres. 
STE 
ERR 
A, L IZRS-O-A 
M ASCO Pog PES Z5 


三 角 测量 法 中 LASCO、STEREO、CME 前 沿 (418). AMA (白色 圆 盘 ) 之 间 的 相对 位 置 ， 以 及 标 
记 这 些 位 置 所 用 到 的 方位 角 和 俯仰 角 " 


结构 或 CME 的 核心 ) 或 自动 使 用 局 部 相关 跟踪 (local correlation tracking, LCT) 方法 来 进 


fr, 一旦 找到 像素 之 间 的 对 应 关系 ， 则 沿 该 像素 视线 方向 反映 到 太阳 ， 最 后 计算 对 应 点 的 三 


维 坐 标 就 可 以 实现 三 维 重 构 。 由 于 视线 必须 在 同一 对 极 面 上 ， 它们 在 这 个 平面 上 的 交点 是 明 


人 确 


的 ， 与 选择 的 目标 无 关 ， 所 有 的 对 极 面 都 在 连接 两 个 航天 器 的 直线 上 相交 。 这 种 方法 通常 


RA see”, BI TP E” 


STEREO B 

Cs 的 < 人、 

pm — Coane Fons 
STEREO A oa 
( 右 眼 图 像 ) 


图 13 ”两 个 航天 器 对 一 个 日 园 特 征 的 三 角 测量 示意 图 


[7 


Srivastava 等 人 外 对 2007 年 5 月 20 日 爆发 的 一 个 部 分 晕 状 CME 进行 了 TP 重 构 。 
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图 14 a) 运动 的 特征 在 COR1B 图 像 上 的 投影 位 置 示意 图 ; b) 包含 运动 特征 的 极 面 示意 图 一 


it STEREO/COR1 的 偏振 图 像 获得 了 总 亮度 图 像 ， 然 后 结合 总 亮度 图 像 和 COR2 的 白光 
图 像 进 行 重建 。 但 在 重建 前 需要 先 对 图 像 进 行 校正 处 理 ， 使 得 STEREO-A 和 STEREO-B 
图 像 的 分 辩 率 以 及 其 中 的 太阳 中 心 坐 标 都 相同 。 这 时 就 可 以 利用 太阳 物理 软件 SSW 中 的 
SECCHI 程序 包 ， 重 建 出 前 沿 上 的 某 一 点 的 三 维 坐标 ， 并 在 三 维 空间 中 确定 前 沿 所 选 点 的 
真实 高 度 和 速度 ， 而 重建 的 真实 高 度 比 投影 高 度 高 1.80 倍 。 

TP 方法 假设 仿 射 几何 在 日 园 仪 中 是 有 效 的 。 这 个 假设 取决 于 太阳 和 目标 之 间 的 距离 与 
太阳 和 观测 者 之 间 的 距离 的 比值 ， 这 个 比值 越 小 ， 假 设 越 合 理 。 通 常 要 重建 的 对 象 与 观测 者 
之 间 的 距离 大 于 200 Ro。， 这 个 距离 比重 建 目 标的 尺寸 和 该 目标 到 太阳 的 距离 要 大 得 多 。 仿 
射 几何 假设 观测 者 位 于 无 限 远 的 位 置 ， 因 此 所 有 的 观测 视线 可 以 被 认为 是 平行 的 ， 与 目标 在 
不 同 观 测 角 度 所 对 应 天 空 平 面 上 的 投影 之 间 的 距离 hh 无关。 这 个 假设 会 由 于 /200 Ro 不 等 
于 零 而 有 误差 ， 这 个 误差 会 导致 重 构 结果 中 的 目标 尺寸 偏 高 。 另 外 ，TP 重 构 方法 本 身 也 存 
在 一 定 的 误差 ， 这 个 误差 取决 于 两 个 航天 器 的 分 离 角 度 ， 分 离 角度 合适 时 ， 误 差 就 比较 小 。 

TP 重 构 衍生 出 很 多 种 重 构 方 法 ， 包 括 可 以 精确 计算 CME 三 维 速度 和 传播 方向 的 三 维 
高 度 时 间 重 构 技 术 (3D-HTS )， 局 部 相关 跟踪 加 上 三 角 测量 的 LOT-TP 技术 ”和 使 用 两 个 
或 三 个 观测 仪器 组 合 (STEREO 或 STEREO + SOHO) 进行 三 角 测 量 的 方法 至 号 。TP 重 
<u FIER] 3D-HT 技术 也 是 建立 在 仿 射 几何 原理 成 立 的 基础 上 ， 根 据 几 何 计算 ( 见 图 14) 获 
= 得 重建 点 在 球 坐 标 系 中 的 坐标 (Rap, 0, 入， 其 中 0 和 入 分 别 是 纬度 和 经 度 。 用 这 种 方法 可 
PT 以 推断 出 CME 在 每 一 个 时 刻 的 位 置 和 传播 方向 ， 其 实际 速度 可 以 通过 估计 Rap 的 时 间 导 
i 数 来 计算 。 

Mierla 等 全 利用 3D-HT 技术 从 两 幅 STEREO 观测 到 的 图 像 中 获得 CME 中 一 个 能 
够 被 很 好 识别 的 特征 的 高 度 - 时 间 图 ， 由 此 得 到 了 两 个 独立 的 投影 速度 矢量 ， 并 可 进一步 构 
造 一 个 三 维 速度 矢量 。HT 技术 的 广泛 应 用 已 经 有 些 时 日 ， 但 主要 用 于 确定 CME 在 天 空 平 
面 的 投影 速度 ， 即 在 日 冕 仪 获 得 的 图 像 中 选择 一 个 特定 的 目标 ， 通 过 跟踪 其 位 置 随时 间 的 变 
化 来 获得 CME 的 传播 方向 。 现 在 由 于 STEREO 双星 数据 的 出 现 ， 使 得 HT 技术 不 再 局 限 
于 求解 CME 的 二 维 速度 。Sheeley 等 人 ” 开发 出 可 以 自动 检测 图 像 中 最 微弱 的 移动 特征 的 
算法 ， 使 得 3D-HT 技术 不 仅 可 以 即时 显示 CME 的 运动 和 速度 ， 还 可 以 计算 模糊 目标 的 三 
维 速度 并 且 灵 敏 度 很 高 。 
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2.1.7 表面 轮廓 拟 和 方法 
Feng 等 人 ”在 2012 年 提出 一 种 基于 多 人 台 日 冕 仪 得 到 的 CME 外 围 反 向 投影 来 重建 


CME 的 新 方法 ， 称 为 表 1 


置 、 形 状 及 


果 更 精确 ， 


并 且 与 了 


轮廓 拟 合 (mask fitting, MF) 方法 ， 其 目的 是 重 构 三 维 CME 的 位 
其 形状 的 演化 。MEF 方法 使 用 的 是 三 个 视点 的 日 冕 仪 数据 ， 比 两 个 视点 的 重建 结 


E 演 建 模 不 同 的 是 MF 方法 不 需要 假设 CME 的 几何 形状 。 


传统 三 角 测 量 法 或 者 点 


t 轿 方法 都 是 在 一 对 日 蝎 仪 观测 图 像 中 选择 一 个 相同 特征 点 ， 


之 后 沿 不 同 航 天 器 与 太阳 之 间 的 连 线 往 太阳 方向 投影 ， 然 后 得 到 它们 的 三 维 坐 标 。 而 MF 方 


向 的 投影 都 


法 是 在 包含 太阳 的 3D 空间 中 选择 一 个 点 分 别 投影 到 三 个 日 园 仪 观测 图 像 中 ， 如 果 沿 视线 方 
立 于 定义 的 CME 区 域 中 ， 则 认为 该 点 是 CME 云 内 的 一 点 ， 再 依次 找到 CME 


边缘 上 所 有 的 点 。 为 了 更 好 地 描述 CME 云 的 三 维 结构 ， 重 建 的 三 维 CME 区 域 必须 足够 接 


近 真 实 CME 
MEF 方法 
RAMEN A 


的 边界 。 
的 具体 步 又” 如下: 首先 在 来 自 STEREO-A, LASCO, STEREO-B 的 三 幅 日 
像 上 分 别 创建 一 个 CME 区 域 ， 将 CME 边缘 曲线 内 的 像素 值 设 为 1， 边缘 曲 


线 外 的 设 为 0;， 然 后 将 一 个 以 太阳 为 中 心 的 三 维 立 方 体 离散 化 ， 其 中 的 每 个 点 分 别 被 投影 到 


STEREO-A, STEREO-B 和 SOHO 获得 的 上 述 三 个 图 像 平 面 上 。 只 有 那些 被 投影 到 三 幅 图 


iy CME 区 域 上 的 点 才 被 认为 是 属于 该 CME 的 点 。 然 后 ， 使 用 Bézier 曲线 平滑 边界 ， 获 
得 如 图 15 所 示 的 CME 的 三 维 边 界 形状 。 


MF 方法 不 需要 假设 CME IAE 
适应 于 各 种 不 规则 的 CME 形状 ， 使 用 起 来 方便 灵活 ”。 基 于 MP 方法 获得 的 CME ZAER 
状 使 我 们 可 以 分 析 其 几何 中 心 和 主攻 


10 
5 
o 
~ 0 
N 
-5 
a 
5 
bp 


TE: 黑 球 代表 太阳 ， 


绿 点 代表 CME 云 的 几何 中 心 。 
图 15 “利用 表面 轮廓 拟 合 方法 重 构 的 三 维 CME z^ 


R， 能 够 很 好 地 重 构 CME 表面 ， 获 得 表面 精细 结构 名 


Ws 


r 


尺度 等 。 该 方法 的 不 足 之 处 在 于 所 得 的 结果 没有 包括 


H 
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CME 的 内 部 结构 ， 同 时 也 要 求 三 个 观测 仪器 之 间 具 有 合适 的 分 离 角 度 。 
2.2 ”基于 日 球 成 像 仪 观测 数据 的 重建 
2.2.1 几何 三 角 测量 法 

前 面 介 绍 的 三 角 测 量 方法 以 及 点 共 斩 方 法 都 可 以 确定 CME 传播 方向 和 径 向 距离 ， 但 这 
两 种 方法 以 及 一 些 其 他 类 似 的 方法 天 ”都 需要 识别 和 跟踪 来 自 两 个 航天 器 的 图 像 对 中 的 相 
同 特征 ， 这 在 距离 较 远 的 地 方 (CME 信号 变 得 非常 微弱 和 弥散 ) 是 不 可 能 的 。Liu S AUT 提 
出 一 种 基于 成 像 观测 重 构 的 时 间 - 延 展 角 图 的 几何 三 角 测 量 法 。 这 个 方法 的 优点 在 于 它 基 于 
时 间 - 延 展 角 图 ， 可 以 首次 将 几何 三 角 剖 分 应 用 于 HI-l1 和 HI-2 的 弱 特 征 ， 就 可 以 在 日 球 层 
中 连续 跟踪 CME 并 预测 其 对 地 球 的 影响 。 在 此 基础 上 ，Liu 等 人 ™ 和 Lugaz 等 人 "将 调和 
平均 值 几何 结构 引入 几何 三 角 测量 框架 中 ，Davies 等 人 "把 自 相似 脱 胀 几何 结构 引入 到 几 
何 三 角 测 量 框架 中 ， 这 些 方法 都 取得 了 成 功 ， 证 明 上 述 几 何 三 角 测 量 的 概念 和 框架 都 是 行 
之 有 效 的 ， 随 后 的 很 多 工作 也 进一步 证 实 了 这 一 点 。 更 重要 的 是 ， 将 这 个 几何 三 角 测 
量 概念 和 实施 方法 与 其 他 方法 进行 比较 ， 如 太阳 风 就 地 测量 、 行 星际 TT 型 射电 暴 频 率 漂移 ， 
所 得 结果 非常 相似 ， 并 进一步 给 出 了 典型 快 、 慢 CME 在 整个 日 地 空间 的 传播 规律 。 

从 原理 上 来 讲 ， 几 何 三 角 测 量 是 基于 图 16， 只 要 沿 着 两 个 航天 器 之 间 的 方向 移动 ， 就 
可 以 看 到 白光 特征 。 特 征 的 延展 角 (特征 相对 太阳 -航天 器 连 线 的 角度 )， 在 STEREO-A 和 
STEREO-B 中 分 别 表示 为 aa 和 we， 可 由 治 黄道 平 面 琶 加 动态 相 减 图 得 到 的 时 间 - 延 展 角 图 
得 到 。 基 于 图 16 简单 的 几何 关系 可 以 得 到 : 


rsin(aa + Ba) 


=d 27 

sin QA EE (27) 

rsin(an + Bp) = (28) 
sin ap 

ba +B =y , (29) 


其 中 , r 是 特征 与 太阳 的 径 向 距离 ，BA 和 Bs 是 特征 相对 于 太阳 -航天 器 线 的 传播 角 ，dA 和 
ds 分 别 是 STEREO-A 和 STRERO-B 与 太阳 的 距离 ，7 是 两 个 航天 器 的 经 度 分 离 角 。 一 旦 
从 成 像 观 测 中 测量 到 延展 角 (wA 和 ap)， 上 述 方程 就 可 以 求解 r+, Ba 和 BBp， 且 解 是 唯一 的 
(与 模型 拟 合 相 比 )。 若 STEREO 航天 器 da ~ dp， 可 得 : 


sin QA sin(arg + y) — sina, sin org 


(30) 


tan 一 一 - , 
PA sin a, cos(arg + y) + cos a, sin ag 


不 过 ，Liu AU" Jah, BME da A dg 也 可 以 快速 获得 CME 的 传播 方向 。 

几何 三 角 测 量 法 相对 可 靠 ， 并 且 没 有 自由 参数 ， 唯 一 的 假设 是 同样 的 特征 可 以 在 两 个 航 
天 器 的 时 间 - 延展 角 图 中 追踪 到 ， 这 样 即使 是 微弱 特征 也 可 以 追踪 到 1 AU 附近 。 这 个 方法 
有 三 个 优点 ”: 首先 ， 基 于 时 间 - 延 展 角 图 ， 首 次 将 几何 三 角 剖 分 应 用 于 HL1 和 HI-2 的 弱 
特征 ; 其次， 与 单轨 道 拟 合 技术 相 比 ， 该 方法 所 依赖 的 假设 量 更 少 ， 因 此 解 的 精度 更 高 ， 第 
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STEREO B 


gi 


图 16 ”几何 三 角 测量 法 基于 的 关系 图 


， 它 可 以 确定 CME 特征 (或 其 他 白光 特征 ) 从 太阳 一 直到 1 AU 的 传播 方向 和 真实 距离 。 


此 方法 的 缺点 是 无 法 获得 CME 三 维 几何 形状 和 磁场 信息 。 


2.2.2 ” 自 相 似 膨胀 方 ; 


由 Davies 等 人 "提出 的 自 相似 膨胀 (self-similar expansion, SSE) 方法 利用 HI 数据 进行 


CME Xf. SSE 方法 是 通过 单 视角 成 像 来 获得 CME 的 时 间 - 延 展 角 轮廓 ， 根 据 获得 的 结果 
探究 CME 的 传播 方向 和 传播 速度 。 早 在 1999 年 ，Sheeley 等 人 ” 提出 一 种 固定 $ AME 


(FPF) 获得 CME 时 间 - 延 展 角 轮 廊 的 方法 。 还 有 一 种 相似 方法 : 调和 平均 值 拟 合 (harmonic 
mean fitting, HMF)” ^, FPF 拟 合用 到 CME 是 一 个 径 向 传播 的 点 源 且 传 播 方向 不 变 的 假 
Y (图 17 a) 的 空心 黑 点 ); HMF 拟 合 则 用 到 CME 的 前 沿 为 贺 形 ， 传 播 过 程 中 膨胀 的 圆 固 定 
在 太阳 中 心 ， 并 且 CME 沿 一 个 固定 的 径 向 轨迹 传播 的 假设 ( 见 图 17 a) 的 实心 黑 点 以 及 灰 


色 大 圆 )。 这 两 种 方法 都 是 SSE 方法 的 极端 情况 。 


SSE 模型 中 CME 的 截面 也 是 圆 形 的 ， 在 与 观测 者 感 兴趣 的 位 置 角 所 对 应 的 平面 上 ， 半 


径 随 着 CME 向 外 传播 增 大 ， 


这 个 圆 也 不 再 固定 在 太阳 上 (如 图 17 b) 所 示 )， 而 是 随 着 CME 


的 传播 向 外 移动 ， 使 CME 对 太阳 中 心 所 张 的 角度 是 个 常数 。 因 此 ， 在 SSE 模型 中 ， 径 向 距 


E Rsspa， 即 CME 的 顶点 a 


可 由 下 式 表示 : 


RssEalt) = 


(距离 太阳 最 远 的 点 ， 见 图 17 b) 中 的 实心 黑 点 ) 与 太阳 的 距离 


do sin[ea(t)|(1 + sin A) 
sinfe,(t)+¢]+sinA  ' 


(31) 


KB. do 是 观测 者 与 太阳 的 距离 ，e。 是 太阳 -观测 者 连 线 和 从 观测 者 出 发 到 CME 圆 形 前 沿 
的 切线 之 间 的 夹 角 ， 入 是 同一 平面 内 CME 的 半 宽 角 (0 € Asse < 90°), o 是 太阳 -观测 者 连 
线 和 CME 传播 方向 之 间 的 夹 角 ， 在 CME 传播 过 程 中 是 被 默认 不 变 的 。 同 理 ，CME 尾 
点 5 ( 即 最 接近 太阳 的 点 ， 图 17 b) 中 的 空心 黑 点 ) 与 太阳 的 距离 Rose 也 可 被 推导 出 ” : 
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if if 
25 观测 者 观测 者 
a) b) 


ik: a) FP 和 HM 模型 的 示意 图 ，FP 用 空心 黑 点 表示 ，HM 用 大 圆圈 /实心 黑 点 表示 ，CME 沿 径 向 方向 传播 ; 
b) SSE 模型 示意 图 ， 空 心 黑 点 表示 CME 的 尾 点 ， 实 心 黑 点 表示 CME 的 顶点 。 


17 ”三 种 方法 模型 示意 图 


do sin[es(t)|(1 — sin A) 
sinfe,(t) + ¢]—sinA ' 
CH, ey 是 太阳 -观测 者 连 线 与 从 观测 者 出 发 的 CME 后 方 切线 之 间 的 夹 角 。 对 于 CME 而 
寺 ， 后 方 与 前 沿 不 同 ， 它 不 太 可 能 被 前 方太 阳 风 的 相互 作用 所 扭曲 ”。 对 这 些 表达 式 的 推 
可 以 参见 文献 加 。 在 上 述 公 式 中 ， 当 入 = 0° 时 ， 得 到 FP 重 构 ， 即 


RssEb lt) = 


(32) 


x 


dm mil 


.. dosin[e(£)] 
Rer(t) = sinfe(t)+ d] ' vn 


m4A=90° 时 ， 我 们 得 到 HM 重 构 ， 即 


2do sin[e(t)] 
sin[e(t) ++ 9] +1 


Rum (t) = (34) 


我 们 注意 到 无 论 是 公式 表述 ， 还 是 模型 的 几何 构图 都 表明 FPE 方法 没有 区 分 CME 的 
前 沿 和 后 方 ， 而 HMF 方法 中 的 CME 后 方 永 远 位 于 日 心 ， 所 以 这 两 种 方法 是 SSE 方法 的 两 
个 极限 情况 。SSE 方法 的 实现 需要 进一步 计算 角度 e (t) Al ey(t)™, Bn 


_, [ —bc + ava? + b? — cP? 
e(t) = cos! ( "NW ) , (35) 
aN +, d 1 ) 
Ttc ; . 
ü dg 6) — cos(Óssg) , b = sin(óssg) , C= csin(Assg) . (36) 
VssEt 


Hop, e 的 符号 按 这 样 的 方式 确定 : ER ealt) 时 ， 取 “二 ”号 ， 而 “一 ”号 时 得 到 e). 
SSE 方法 通过 单 视角 观测 到 的 CME 瞬 态 的 时 间 - 延 展 角 轮廓 拟 合 到 这 个 方程 ， 以 获得 径 向 
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速度 和 传播 方向 的 最 佳 拟 合 。 从 这 些 参数 可 以 导出 我 们 感 兴趣 的 其 他 参数 ， 如 预测 到 达 指 定 
位 置 的 时 间 。Davies 5& A" 使 用 蒙特 卡 洛 (Monte Carlo) 模拟 来 验证 SSE 模型 的 适用 条 件 ， 
得 到 延展 角 应 在 19° ~ 74° 之 间 ， 过 大 或 者 过 小 的 延展 角 都 会 导致 结果 不 可 信 。 而 延展 角 不 
能 过 大 则 将 限制 SSE 模型 对 空间 天 气 预 报 的 适用 范围 。 由 于 SSE 模型 假设 CME 均 有 一 个 
圆 形 的 瞬 态 截面 ， 明 显 与 实际 情况 有 些 出 入 ， 重 构 的 结果 与 真实 的 CME 位 形 肯 定 有 差距 。 
目前 研究 人 员 正 在 努力 优化 该 技术 ， 使 之 适用 于 具有 椭圆 形 瞬 态 截面 的 CME 重 构 (如 2.2.3 
节 介绍 的 椭圆 演化 模型 )。 
2.2.3 ”椭圆 演化 模型 
2014 年 1 月 7 日 ,一 个 面向 地 球 的 活动 区 爆发 了 一 个 快速 CME (天 空 平面 的 投影 速度 
约 2400 km/s)。 由 于 从 靠近 太阳 中 心 的 源 区 爆发 的 快速 CME 通常 会 影响 地 球 ”™， 所 以 
许多 观测 者 预测 这 次 CME 会 对 地 球 产 生 很 大 的 影响 ， 然 而 ， 随 后 并 没有 发 生地 磁 暴 。 这 说 
明 此 次 事件 并 没有 沿 径 向 传播 ， 并 且 这 种 非 径 向 的 传播 方向 是 在 很 接近 太阳 的 地 方 就 确定 
了 ， 而 不 是 在 行星 际 空间 传播 过 程 中 发 生 偏转 所 致 ”。 为 了 研究 此 次 CME 这 样 传播 的 具体 
Y ER, Móstl A 提出 一 种 新 的 椭圆 演化 (ELlipse Evolution, ELEvo) 模型 ， 专 门 针对 此 
> 次 快速 CME 的 非 径 向 传播 特征 。 

t od ui 
型 中 的 解析 表达 式 可 以 将 激 波 的 形状 投影 到 给 定 的 行星 或 局 地 观测 航天 器 的 天 空 平 面 上 ， 
ps i uc x mu P 从 
日 冕 仪 观测 得 到 的 CME 初始 速度 、 方 向 和 宽度 可 以 作为 初始 条 件 。 


CME 中 心 方向 


ik: a) CME 的 前 沿 、 激 波 或 磁 通 量 绳 的 前 端 ， 被 描述 为 一 个 椭圆 ， 在 黄道 面 或 太阳 赤道 面 以 恒定 的 角 宽 和 长 短 
轴 比 向 太阳 外 传播 ，b) 推导 CME 边界 任意 点 的 速度 的 几何 结构 。 


18 ”ElEvo 模型 中 自 相 似 扩展 椭圆 的 几何 推导 


ELEvo 模型 用 到 如 下 假设 : (1) CME 边界 的 日 球 经 度 角 宽度 一 直 保 持 不 变 ，(2) 椭圆 
( 激 波 ) 的 一 条 主轴 沿 传播 方向 ，(3) 椭圆 的 长 轴 和 短 轴 之 比 为 常数 ，(4) 椭圆 对 称 轴 的 方向 
不 变 。 对 于 在 行星 际 传播 的 CME 则 采用 DBM 模型 ” 。 所 以 包含 DBM 的 ELEvo 模型 涉 
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及 到 4 个 自由 参数 ， 它 们 分 别 是 : (1) 长 短 轴 之 比 的 倒数 fs (2) 椭圆 的 半角 宽度 A; (3) BA 
力 参 数 7 (AA (0.1 ~ 2)x107 km ^P): (4) 背景 太阳 风速 度 w。 使 用 DBM 模型 时 选择 较 小 
的 阻力 参数 可 以 描述 距离 太阳 在 1 AU 以 内 的 激 波 传播 ， 从 而 可 以 计算 CME 激 波 到 达 地 球 
的 时 间 和 速度 。 

我 们 现在 推导 ELEvo 模型 中 所 用 到 的 公式 以 及 沿 椭圆 前 沿 的 每 个 点 的 速度 只 。 椭 
顶点 ( 沿 椭圆 中 心 方向 离 太 阳 最 远 的 椭圆 点 ) 的 R(t) FH DBM™ 给 出 ， 长 短 轴 之 比 的 倒数 
f —b/a, Ti a, — Qa/b， 可 以 得 到 如 下 结果 : 


Ted 


f =b/a 
B-—A 3 (37) 
0 = tan ! (b? /a? tan f) 


其 中 ， 角 5 是 半 长 轴 a 5A T WIEDER ZA ( 见 图 18 a)), T 的 位 置 是 一 条 以 太阳 
为 起 点 的 直线 在 椭圆 上 的 切 点 ， 由 几何 关系 可 以 看 出 8 = A. 0 是 椭圆 的 极 角 ，8 和 0 之 间 
的 几何 关系 由 式 (37) 给 出 ，7 是 椭圆 中 心 与 工 点 的 距离 ， 结 合式 (37) 得 到 : 


0 —tan !(f^tanf) . (38) 


从 图 18 a) 中 橙色 大 三 角形 上 的 正弦 定律 ， 推 导出 : 


sin 入 sin @ 


r R@b’ (39) 
其 中 ， 角 a 由 三 角形 的 几何 关系 得 : a = 90° 十 9 一 入 则 距离 7 在 极 坐 标 系 下 的 定义 为 : 
r = b/V/(f2—1)cos?O+1 . (40) 
将 a 和 7 的 方程 代入 式 (39)， 并 令 w = (fF? — 1) cos? 6 +1 可 得 椭圆 的 a,b,c 参数 ; 
= a=6/f 
~ = e : ^ J on 
c= R(t) — b 


式 (41) 为 ELEvo 椭圆 模型 参数 的 最 终 描述 。 从 式 (38) 和 (40) 可 知 椭圆 的 短 轴 b 依赖 
于 所 有 已 知 变量 IG, f, Ap fia 则 简单 地 遵循 式 (37) 中 f 的 定义 。 得 到 了 模型 的 几何 
状 之 后 ， 就 可 以 计算 椭圆 上 任意 点 的 速度 。 为 了 计算 方便 ，Mistl 等 人 ™ 在 椭圆 中 引入 一 
个 笛 卡 尔 坐标 系 (图 18 b)， 其 中 X HEA CME 传播 方向 ，Y EXT X, c 是 从 太阳 中 
心 到 椭圆 中 心 的 向 量 ，d 是 一 条 沿 太阳 -地 球 连 线 的 直线 ， 终 点 在 椭圆 边界 (从 图 18 b) 可 以 
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看 出 4 与 椭圆 有 两 个 交点 )，" 连接 椭圆 中 心 和 d 在 椭圆 上 的 终点 ，A 是 CME 中 心 方向 和 太 
阳 - 地 球 线 之 间 的 夹 角 ， 地 球 也 可 能 是 其 它 行星 或 者 航天 器 。 根 据 上 述 公 式 ， 可 以 计算 出 
与 椭圆 的 两 个 交点 : 


ccos A + "C — c?) f? sin? A + b? cos? A 
di2 = 2.: 2 E (42) 

f? sin^ A+ cos? A 
其 中 ， 根 号 前 带 正 号 的 解 是 “前 沿 ” 解 (di1)， 负 号 的 解 是 “后 方 ” 解 (d2)。 椭 圆 在 角度 所 定 
义 的 位 置 上 的 速度 是 由 椭圆 的 自 相 似 膨 胀 推导 出 来 的 ， 这 意味 着 半 宽 不 变 。 自 相似 膨胀 的 
假设 意味 着 CME 的 形状 不 会 随时 间 改 变 ， 因 此 ， 椭 圆 上 所 有 点 的 速度 与 距离 之 比 必须 为 常 


数 ”: 
_ a(t) 
Va(t) = Ri V(t). (43) 
此 外 ， 当 dy (t) 等 于 行星 的 日 心 距离 时 ， 局 地 观测 者 即 可 以 给 出 椭圆 到 达 某 个 局 地 位 置 
T 的 时 间 ， 以 及 相应 的 速度 Va(t). HF ELEvo 模型 中 椭圆 的 长 短 轴 之 比 是 一 个 自由 参数 ， 所 
^ 以 对 椭圆 形状 的 控制 比 自 相 似 膨 胀 圆周 更 灵活 ， 更 适合 多 点 局 地 观测 的 自 兴建 模 。M6stl 等 


人 四 利用 地 球 和 火星 多 点 局 地 观测 的 数据 限制 了 ELEvo 模型 的 形状 (长 短 轴 比 为 1.4 土 0.4)， 
KING Bl KE CME 运动 的 方向 ”。 这 些 结果 证 明 这 次 的 CME 事件 没有 引起 地 磁 暴 
不 是 因为 CME 在 行星 际 空间 的 偏转 导致 ， 而 是 由 于 CME 源 区 附近 的 强 磁场 导致 CME 一 
开始 的 非 径 向 传播 。 由 于 观测 到 的 CME 图 像 都 是 沿 视线 积分 的 结果 ， 所 以 CME 在 日 球 经 
度 上 的 非 径 向 运动 的 研究 很 难 进行 ， 而 ELEvo 模型 则 为 研究 CME 的 非 径 向 传播 提供 了 一 
种 有 效 的 方法 。 
2.24 相关 辅助 重 构 法 

基于 日 球 成 像 数 据 的 重 构 方法 还 有 一 种 相关 辅助 重 构 (CORrelation-Aided Reconstruc- 
tion, CORAR) 法 ”， 这 个 方法 可 以 自动 识别 和 定位 太阳 风 的 非 均 匀 结 构 ， 特 别 是 小 尺度 瞬 
变 结 构 。 后 来 ，Li 等 人 改进 了 这 个 方法 ， 使 得 该 方法 可 以 应 用 到 三 维 OME 重 构 ， 并 从 
STEREO/HI 的 双 图 像 中 自动 识别 沿 日 地 连 线 传播 的 CME。 与 其 他 方法 相 比 ，CORAR 77 
法 简单 很 多 ， 并 且 对 所 感 兴趣 的 目标 没有 任何 特殊 的 形态 假设 。 

CORAR 方法 首先 需要 对 HI 图 像 进行 处 理 ， 以 消除 FF ANB RB BNO. FAN 
天 器 位 置 、 视 场 和 基线 (日 地 连 线 ) 的 情况 下 ， 根 据 两 个 STEREO 航天 器 在 日 心 -地 球 -黄道 
(HEE) 坐标 系 下 的 时 间 空 间 位 置 ， 将 HI-1 图 像 投 影 到 基线 所 在 的 子午 面 上 。 然 后 用 纬度 和 
到 太阳 中 心 的 距离 作为 极 坐标 ， 对 图 像 沿 径 向 投影 ， 并 在 这 个 子午 平面 上 标 出 CME 的 空间 
位 置 。 使 用 一 个 径 向 纬度 抽样 框 来 获得 一 个 二 维 局 部 亮度 变化 ， 然 后 沿 “ 基 线 ” 计 算 相 关系 
数 cc。 如 果 目 标 在 基线 上 ， 这 两 个 部 分 是 很 好 的 重 奴 ( 见 图 19 a))， 那 么 它们 应 该 是 高 度 相 
KN; WR A eh AAR, SARA SR) et LEAS ( 见 图 19 b))， 则 相 
关 性 较 低 。 

下 面 介绍 如 何 将 HI 图 像 投影 到 子午 平面 上 。 首 先 假设 P 是 经 度 为 o 的 子午 平面 上 像素 
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AIH e Ane 
hw Pd Seal cd He 
Dy | ee 
ZS pope s 
P ad "i E dt STEREO-A 


STEREO-B a l 
a) b) 


VE: MEA BA EAR 2 nee EER. REINER. Æ STEREO-B/STEREO-A 上 的 两 条 细 绿 / 蓝 
径 向 实 线 是 恰好 与 该 物体 相 切 的 视线 ， 而 粗 绿 / 蓝 实 线 是 覆盖 HI-1 图 像 特征 物体 的 基线 的 分 割 。a) 当 特 征 对 象 
基线 上 的 情况 ，b) 当 这 个 特征 对 象 远离 基线 的 情况 。 


19 ”CORAR 法 原理 图 名 


的 投影 点 的 HEE 坐标 ， 则 有 : 


P(z,y,z)=O+kn ， (44) 


KH, n 是 HI-1 图 像 中 像素 相对 于 观测 者 O 的 方向 ， 可 以 从 STEREO/HL1 的 头 文 件 中 获 
HU 大 是 未 知 的 正 数 。 投 影 点 的 纬度 入 和 与 太阳 中 心 之 间 的 距离 D 是 投影 过 程 中 的 两 个 关 
键 未 知 参 数 ， 它 们 之 间 的 关系 如 下 : 


| x = DcosAcosó 


y = DeosXsing , (45) 
z= DsinA 


其 中 ，% 是 投影 点 的 经 度 。 通 过 求解 上 述 两 个 方程 可 以 计算 出 xz,y,z,k, 入 和 DD。 根据 导出 的 
AMD, EZEN o 的 子午 面 上 绘制 两 幅 HI 的 投影 像 。 在 同一 子午 面 上 ， 两 幅 投 影 图 像 的 
重 登 区 域 的 相关 系数 cc。 先 划 出 一 个 大 小 合适 的 动态 3D ( 径 向 -纬度 -时 间 ) 区 域 (纬度 5°, 
径 向 约 8 Ro); ee a p di 用 于 cc 的 计算 。 采 样 区 域内 的 亮度 变化 主要 取 
RF CME 密度 的 空间 变化 。 为 了 突出 密度 空间 变化 规律 ， 采 样 区 域 的 亮度 数据 减 去 它们 的 
平均 值 。 相 关系 数 cc 的 计算 是 由 线 " Pearson 相关 系数 给 出 : 


Da 1 — p)(ai = q) (46) 
^ Ju 9n 


其 中 ，p 和 g 两 组 数据 集 (对 应 于 数据 在 同一 取样 箱 )， 样 本 大 小 为 mw, 石和 7 分 别 是 和 的 
平均 值 。 在 我 们 选取 的 抽样 区 域 中 n= 11 x 41 x t,t 等 于 时 间 步 长 的 数量 ， 可 以 是 1 也 可 
以 是 5. cc 的 取 值 范围 为 -1 ~ 1: cc 值 越 大 ， 相 关 性 越 好 ; cc 值 为 负 则 认为 是 不 相关 的 。 

X (46) 表明 ， 只 要 信 品 比 足够 大 ，cc 的 值 不 受 特征 密度 的 影响 ， max 密度 的 变化 规 
售 的 影响 CME 的 位 置 、 信 噪 比 等 因素 可 能 会 影响 计算 出 的 cc 值 ， 需 要 进一步 研究 和 修 
IE . 20 展示 了 CORAR 方法 对 一 个 CME 重 构 的 结果 。 为 了 验证 此 方法 的 可 靠 性 ，Li 
等 全 将 GCS 方法 应 用 到 同一 个 CME 事件 ， 结 果 发 现 这 两 个 方法 得 到 的 CME 的 角 宽 度 
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和 传播 方向 基本 一 致 ， 不 同 之 处 则 是 在 CME 的 几何 结构 上 ， 从 图 
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方法 得 到 的 CME 的 结构 不 再 是 传统 的 磁 通 量 绳 形状 ， 更 像 是 
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20 可 以 看 出 由 CORAR. 
个 非常 扭曲 的 结构 。 


2010-04-03 16:49:01 


地 球 


注 : 最 左边 显示 的 是 三 维 采 样 盒 中 的 高 cc 区 域 。 从 左 到 右 是 从 不 同 
JE y 轴 方 向 在 )， 地 球 ( 沿 着 x 轴 方 向 在 HEE 坐标 ) 和 顶部 ( 沿 着 z 


20 用 CORAR 方法 对 一 个 CME 


CORAR 方法 的 优点 是 不 需要 假设 CME 几何 形态 ， 


通过 Li E AC 的 改进 ， 此 方法 还 可 以 得 到 三 
法 也 存在 一 些 不 足 。 如 共 面 效应 ， 它 是 指 如 果 CME Mg 


部 分 ( 见 图 
得 到 相似 的 投影 


图 像 ， tcc 值 较 高 ， 这 会 导致 “ 


EX] JH 


三 维 空间 中 CME 的 径 向 速 


度 非常 低 ， 而 大 尺度 CME 的 重建 会 偏向 于 两 个 航天 器 之 间 的 中 
追踪 位 于 预 设 基 线 附近 的 CME， 而 不 是 整个 HI 视 场 中 的 CME。 此 外 ，CORAR 方法 和 三 


mim 


Tii 


o 
FE 
m e 
i, à 


度 看 到 的 CME， 观 测 方 向 分 别 是 东 ( 沿 着 
册 方 向 )。 黄 色 的 点 是 太阳 。 


可 以 自动 地 对 CME 进行 重 构 。 并 


度 分 布 。 但 是 此 方 
F 两 航天 器 连接 线 组 成 的 平面 的 中 心 
19 a))， 即 两 航天 器 视线 相同 ， 那 么 即使 CME 远离 基线 ( 见 图 19 b))， 我 们 
k 面 ”附近 CME 的 经 
FP 心平 面 。CORAR 方法 只 能 


也 会 


度 方向 的 定位 精 


角 测 量 方法 一 样 ， 都 使 用 简单 的 光学 投影 ， 很 可 能 将 大 尺度 CME 的 两 个 不 同 部 分 错误 地 视 


为 同一 部 分 所 。 
2.3 小结 


因此 ，CORAR 方法 对 于 小 上 


度 瞬 变 定位 更 准确 。 


章节 描述 了 10 种 不 同 的 CME 三 维 重 构 方法 ， 在 此 对 它们 使 月 


得 到 的 参数 以 及 方法 的 特点 总 结 在 表 1 当中 。 


3 CME 驱动 激 波 重 构 方 法 


HEIZ 


与 CME 成 像 重 构 具有 较 高 关联 度 的 是 对 CME 驱动 的 激 波 的 三 维 重 构 ， 也 是 近 几 年 来 


新 的 重要 发 


REA. CME 爆发 后 能 量 不 仅 可 以 以 加 速 粒 子 的 
情况 下 ， 快 速 磁 声 波 被 观测 为 Moreton 波 和 EUV XX 


II 型 射电 暴 ””。 当 CME 速度 相 比 于 周围 太 F 
度 时 就 会 在 CME 的 前 方 产生 激 波 。 激 波 在 昌 


风速 度 


图 像 吕 


P3 
In] 


[X 


Nea 


Fa 


沿 形态 ”。 激 波 重 构 的 重要 性 在 于 ， 本 
RETT 、 高 能 粒子 特性 ”"“、 太 阳 风 就 地 观测 


api 


式 释 放 ， 还 可 以 以 波 的 


EX 
EDT I 和 


于 当地 快 磁 声 波 速 或 者 Alfvén XE 
常 表现 为 弓形 激 波形 态 或 者 双 前 
以 给 出 激 波 的 三 维 结构 和 运动 学 特征 ， 并 可 与 CME 


ER, a ae SE PES 


间 天 气 预报 和 高 


ChinaXiv& (ERAT 


| 


> p f 


"To ^P 
L| 2 U c w 


IV: 


Va 


AN 


ch 


180 天 文学 进展 40 卷 
表 1 CME 三 维 重 构 方法 对 比 总 结 
方法 ” ”使 用 数据 ” ”所 需 视角 反 演 方式 得 到 参数 特点 
AI 日 黑 仪 数据 —— 单 视角 正 演 ETE EHER CME 
等 效 散 射 中 心 与 天 空 ee ee sp a 
ARIE BANAWE 单 视角 。 反 演 平面 的 真实 距离 以 及 IN 
CME 的 三 维 形 状 eee ee á 
ET = see, CME 源 区 位 置 、 前 沿 ”对 单个 CME 重 构 的 三 维 速度 
GOS RUA WAA ER 三 维 速度 和 传播 方向 “和 传播 方向 的 估计 非常 准确 
园 仪 或 日 球 。 sus CME 在 1 AU 处 的 。 ”从 结果 中 可 以 推断 出 磁 通 
3DCORE 成 像 仪 数 据 FRA ER 速度 和 磁场 信息 量 强 中 磁场 的 南 向 分 量 
Seti CEN " CME 与 太阳 的 真实 可 快速 估计 CME 源 区 位 置 、 
角 测量 ARDC IARA 。 色 演 距离 与 传播 方向 传播 方向 以 及 修正 投影 效应 
neca 对 于 使 用 局 部 相关 跟踪 技术 
TP BWSR 。 双 视 角 — ER El HEIR 时 应 控制 两 个 视点 的 分 离 
角度 (应 在 30°~150° 之 间 ) 
DM 在 不 假设 CME 结构 下 可 以 
MF RAE 三 视角 mk OANE 很 好 地 重 构 CME 表面 且 三 
状 及 形状 的 演 个 视角 的 重 构 结果 更 精 玉 
m ER T nk EJ ul Je UC I | F 参数 ， 解 更 精确 , 
SSE 求 成 像 仪 数 据 单 视角 。 正 演 。 ” 径 向 速度 和 传播 角度 。 ”延展 角 应 在 19°~74° 之 间 
mem i o 沿 李 圆 任意 点 看 的 速度 “可 研究 快速 CME 的 非 径 向 
ELEvo ”日 球 成 像 仪 数据 单 视角 ER 及 到 达观 测 者 的 时 间 ”传播 特征 
_ CME 的 三 维 结构 及 三 维 可 自动 识别 和 定位 沿 日 地 
CORAR 日 球 成 像 仪 数据 双 视角 SA 空间 中 径 向 速度 的 分 布 ， 连 线 传播 的 CME 


能 粒子 加 速 的 研究 具有 重要 意义 。 此 章 将 简单 介绍 CME 驱动 激 波 的 三 维 重 构 方法 。 
3.1 ” 激 波 球体 模型 


激 波 球状 拟 合 最 早 是 1 


Vourlidas 和 Ontiveros” 提 出， 后 来 Hess 和 Zhang ”也 应 用 


此 方法 来 研究 激流 的 演化 情况 。2017 年 7 月 23 日 爆发 了 一 个 极 高 速 CME (前 沿 速 度 约 
3000 km/s)， 这 个 CME 在 三 个 视点 (STEREO-A, STEREO-B #1 SOHO) 的 观测 中 都 表现 


出 球形 结构 ， 且 速度 极 快 。 


胀 ， 假 设 了 一 个 球 


波 ， 分 别 是 传 


cm 


AP! ge ug 


]TH 


f En, rz 
Liu 等 人 fed 


ix 3 WEE, 


的 三 角 测 量 法 。 如 


Liu 等 人 四 为 了 研究 这 个 晕 状 CME 相关 激 波 的 结构 、 传 播 和 膨 
结构 来 模拟 激 波 (如 图 21 所 示 )。 

假设 这 个 中 心 会 移动 的 球体 激 波 的 传播 方向 不 变 ， 模 型 用 4 个 自由 参数 来 
BOIS FURIA RE, RERE (r) 以 及 球体 中 心 与 太阳 的 距离 (d)。Liu 等 
园 仪 和 日 球 成 像 观 测 的 激 波 传播 和 膨胀 的 距离 计算 方法 ， 此 方法 类 似 于 


述 这 个 激 


图 21 所 示 ， 由 STEREO-A, STEREO-B 以 及 SOHO 


从 三 个 不 同 视角 的 观测 来 确定 距离 的 关系 图 : 
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X/AU 


10L SOHO) 


| 1 I 1 L | 1 1 


1 1 
-1.0 -0.5 0.0 0.5 1.0 
Y/AU 


Hs 灰色 区 域 分 别 是 HI1 在 STEREO-A 和 STEREO-B 上 的 视 场 。 
图 21 ” 激 波 球体 模型 在 黄道 面 的 传播 和 膨胀 


r 


- +d=da " (47) 
sin QA 


d? +r? — 2dr sin(y — ag) = dà sin? ag + d? cos?(y — ap) , (48) 


Et, da 和 dg 分 别 是 STEREO-A 和 STEREO-B 与 太阳 的 距离 ，awA 和 og 则 是 由 
STEREO-A 和 STEREO-B 沿 其 切线 测量 的 激 波 延展 角 ，Y 是 STEREO-A 和 STEREO- 
B 的 经 度 分 离 角 。 激 波 前 沿 与 太阳 的 距离 定义 为 > 十 d， 再 根据 距离 -时 间 图 计算 传播 速度 ， 
从 激 波 球体 的 传播 速度 中 可 以 分 离 出 膨胀 速度 。 图 像 正 向 建 模 给 出 的 距离 和 速度 与 SOHO 
和 STEREO-B 三 角 剖 分 的 距离 和 速度 基本 一 致 ， 但 略 低 于 后 者 ”。 图 22 显示 了 利用 激 波 
球体 模型 分 别 在 三 个 视点 对 CME 激 波 的 重 构 结 果 。 结 果 表 明 ， 无 论 是 日 冕 仪 成 像 重 构 结 果 
还 是 射电 ID 型 暴 和 局 地 观测 结果 ， 此 次 激 波 都 可 以 用 球形 结构 合理 地 模拟 。 


b^ COR? | 


图 22 MORTAR DE S MLAB ama EP 
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3.2 CRT! 

Kwon 等 人 四 为 了 确定 CME 多 个 前 沿 的 三 维 结 
构 和 运动 学 性 质 ， 开 发 了 一 个 复合 模型 。 在 这 个 复 
合 模 型 中 ， 气 泡 状 的 激 波 前 沿用 椭 球 模型 描述 ， 而 
内 部 的 CME 前 沿 则 用 GCS 模型 拟 合 。 为 了 构建 这 
个 复合 模型 ， 定 义 了 两 个 坐标 系 LA 23): 参考 坐 
标 系 (zrerf, Yref, Zret) 和 局 部 坐标 系 (x,y,z). BBA 
标 系 的 原点 Oe MEAP ÙO, zer 轴 是 太阳 的 自 
Erih PAR, Cret 与 地 球 处 看 到 的 太阳 中 央 子 午 线 相 
交 。 在 参考 坐标 系 中 定义 一 个 半径 为 h 的 球体 ， 而 局 
部 坐标 系 的 原点 O 则 位 于 这 个 球体 的 表面 ，z 轴 与 


Em SEER RABE RATE 通 过 OA zne 轴 的 纵向 大 圆 相 切 ，z 轴 朝 着 球体 的 径 


图 四 


Ji x 点 Ox (Tief, E? Wief ,E? Zret, E) 


向 方向 ， 局 域 坐 标 系 与 CME 的 径 向 运动 同步 。 
图 24 显示 了 在 上 面 定义 的 坐标 系 中 创建 椭 球 模 


型 和 GCS 模型 。 椭 球 模型 需要 用 到 7 个 几何 参数 : 3 个 参数 在 参考 坐标 系 中 局 部 坐标 系 的 


或 者 Og(hg, On, dn)» KF he, 0g 和 bp 分 别 是 原点 OF TES 


考 坐 标 系 当 中 的 高 度 、 纬 度 和 经 度 ， 其 他 3 个 参数 是 椭 球 模型 的 三 个 半 主 轴 的 长 度 a, 5 和 


c， 最 后 一 个 参数 ys. per 


如 下 : 


a 与 轴 2 的 夹 角 。 如 图 24 a) 所 示 ， 在 局 部 坐标 系 下 确定 椭 球 


Lp = ACOS[LCOSV, Yp = bcospsinv, z_=csinp , (49) 


其 中 ， 纬 度 0 <p < 180° MAR O< v < 360°. Ha=b= ch], WERKER. GCS 模 
型 与 Thernisien 4 A" 所 用 方法 不 同 之 处 在 于 原点 没有 位 于 太阳 中 心 ， 且 骨架 位 于 r-z F 


面 上 。 


图 24 


a) 激 波 椭 球 模型 ，b) CCS 模型 示意 图 加 
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复合 模型 的 构建 是 一 个 多 次 重复 的 过 程 。 首 先 要 对 自由 参数 进行 初步 猜测 ， 在 参考 坐标 
系 下 建立 椭 球 模型 和 GCS 模型 ， 其 次 ， 计 算 该 结构 在 二 维 平 面 上 的 视图 ， 与 航天 器 观测 至 
的 实际 图 像 作对 比 ， 然 后 将 计算 的 观测 结果 与 实际 观测 到 的 前 沿 进行 比较 ， 最 后 将 参考 坐标 
系 中 构建 的 结构 的 几何 图 形 转换 为 观测 坐标 系 (Xscy Yeo Zee)» TAL ESI Ji A Ose 位 于 
太阳 中 心 ， 轴 Xs。 朝向 观测 者 ， 轴 Yee 和 Zs 在 图 像 平 面 上 分 别 指向 正 西 和 正 北 。 
为 了 重 构 出 三 维 空间 中 的 CME， 必 须 考虑 每 个 像素 点 的 实际 光路 ， 需 要 将 构建 的 三 维 
结构 的 坐标 分 别 转换 为 图 像 平面 向 西 和 向 北 到 图 像 中 心 点 的 角 距 离 (w,v)。 点 (X, Y", 2 ) 在 
观测 坐标 系 中 的 角 距 离 定义 如 下 : 


u = tan”! - y v= tan"! e ; (50) 
Rr, Y"-2Y'—-AY,Z"'—-Z'-AZ,d-D-X'. BRAY M AZ 为 图 像 平 面 上 图 像 中 
心 相对 于 太阳 中 心 的 偏 移 量 。 另 外 ，D 是 观测 者 到 YZ 平面 的 距离 ， 即 


> D = d$ - (AY? € AZ?) , (51) 
CD 

CO Rh. do 是 观测 者 到 太阳 中 心 的 距离 。 反 复 地 重复 这 些 过 程 ， 直 到 重 构 出 的 几何 结构 能 

e 很 好 地 重 现 所 有 仪器 的 观测 结果 。 

= 复合 几何 模型 (包括 椭 球 模型 和 GCS 模型 ) 应 用 三 个 不 同 视角 的 观测 数据 ， 可 以 很 好 地 


再 现 CME 及 其 驱动 激 波 的 三 维 结构 。 对 于 CME 的 三 分 量 结构 ， 最 外 层 微 弱 的 激 波 边界 可 
以 用 椭 球 模型 很 好 地 再 现 ， 而 明亮 的 CME 前 沿 则 用 GCS 模型 再 现 。 两 个 前 沿 与 两 个 不 同 
的 几何 形状 同时 拟 合 表明 ， 它 们 不 是 在 天 空 平 面 上 单个 三 维 结构 的 投影 ， 而 是 具有 本 质 上 不 
- 同 的 三 维 形态 ”。Liu s AUT 利用 激 波 椭 球 模型 重 构 了 CME 和 激 波 的 三 维 几何 形状 ， 并 将 
N ENLIL MHD 模拟 与 局 地 测量 结果 进行 了 比较 ， 发 现 激 波 膨胀 在 日 球 层 激 波 的 范围 中 起 着 
2 重要 的 作用 。 基 于 Liu AT 的 激 波 几何 拟 合 ，Zhu 等 人 "结合 背景 太阳 风 的 稳 态 MHD 
数据 ， 研 究 了 高 能 粒子 特性 ， 并 证 明了 激 波 侧翼 具有 有 效 的 粒子 加 速度 作用 。 

由 于 椭 球 模型 对 激 波 的 形状 有 所 限制 ， 所 依据 此 反 演 得 到 的 结果 难以 给 出 激 波 的 真实 
状 和 结构 。Feng ZAT 结合 对 极 几何 方法 下 ， 对 2.1.7 节 当 中 介绍 的 表面 拟 合 方法 进行 
了 改进 ， 在 研究 分 析 一 个 具体 事件 中 出 现 的 CME 激 波 结构 时 ， 重 构 了 CME 激 波 的 三 维 表 
面 ， 得 到 了 该 表面 形状 随时 间 的 演化 。 他 们 用 激 波 球体 模型 、 椭 球体 模型 以 及 MF 方式 重 
构 了 激 波 三 维 表面 。 三 种 方法 给 出 了 一 致 的 激 波 传播 方向 ， 而 且 MF 方法 还 可 以 捕捉 到 激 
波 波 前 的 外面 结构 。 

Feng 等 人 ”进一步 指出 ， 准 确 重 构 激 波 面 的 三 维 结构 并 得 到 相应 的 激 波 参数 ， 对 理解 
激 波 结构 和 这 些 参 数 与 激 波 加 速 的 太阳 高 能 粒子 之 间 的 联系 至 关 重 要 。 此 外 ， 基 于 不 同 观测 
设备 获得 的 资料 ， 利 用 改进 的 MF 方法 还 可 以 对 我 们 感 兴趣 的 事件 进行 多 方位 的 研究 。 
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4 不 同 重建 方法 应 用 结果 的 比较 


| 由 


在 前 面 的 章节 中 ， 我 们 介绍 了 目前 国际 同行 常用 的 CME 空间 结构 三 维 重 构 的 技术 和 方 
法 ， 讨 论 了 各 自 的 优势 和 仍然 存在 的 可 以 进一步 改进 的 空间 。 我 们 发 现 ， 每 个 重 构 方法 在 一 
定 程 度 上 能 够 较 好 地 重 现 某 一 类 CME 的 某 些 特征 ， 但 都 会 在 其 他 的 方向 与 实际 情况 存在 差 
异 。 受 观测 条 件 和 观测 位 置 的 限制 ， 这 种 差异 还 将 在 以 后 存在 一 段 时 间 。 我 们 在 本 章 中 对 这 
些 重 构 方法 对 具体 事件 的 应 用 结果 做 对 比 ， 对 其 优势 以 及 不 足 做 综合 评价 。 

4.1 2007 年 5 月 15 日 事件 

我 们 讨论 的 第 一 个 事件 是 发 生 于 2007 年 5 月 15 日 的 一 个 CME， 它 的 源 区 AR10956 位 
置 是 N02*E47*，LASCO 给 出 的 CME 前 沿 的 投影 速度 为 491 km/s。 此 时 STEREO 的 两 个 
航天 器 的 分 离 角 度 约 为 8.63°。 

Mierla 等 人 中 利用 3D-HT 方法 得 到 的 CME 速度 约 为 169 km/s, AYN E70°, Bike 
估计 在 N14 左右 。Temmer 等 人 四 利用 TP 方法 ， 通 过 分 析 STEREO-A 和 STEREO-B 
= 的 数据 分 别 获 得 这 次 事件 的 经 度 为 E46。 和 E509, 4 EAP HAE N02 和 NO1*, CME 的 前 沿 
Y 速度 约 为 445 km/s。 这 两 种 方法 得 到 的 结果 差异 较 大 ， 主 要 原因 是 Mierla A 重 构 的 是 
CME 内 部 ， 而 Temmer 4& A" 重 构 的 是 CME 的 边缘 部 分 。 通 过 这 两 个 结果 的 比较 可 以 看 
出 ， 重 构 CME 的 边缘 部 分 更 能 反映 CME 的 综合 信息 和 大 尺度 传播 特征 。 
Mierla 等 人 ”用 另外 三 种 方法 再 一 次 分 析 了 该 事件 ， 推 导 得 出 了 CME 的 传播 方向 和 
特征 的 位 置 (如 表 2 所 示 )。 第 一 种 方法 是 优化 圆柱 壳 方法 (GCS)。 第 二 种 方法 是 先 利 用 局 
部 相关 跟踪 法 (LCT) 确定 STEREO-A 和 STEREO-B 两 张 日 冕 仪 图 像 的 相同 特征 ， 之 后 利 
] TP 方法 重建 CME， 这 种 方法 称 为 LCT-TP 方法 。 第 三 种 方法 是 通过 TP 重 构 求 其 质心 
(CM-TP 方法 ) 的 三 维 坐标 ， 先 假设 汤姆 孙 散 射 是 各 向 同性 ， 则 像素 的 亮度 与 沿 各 自视 线 的 
电子 质量 密度 的 积分 成 正比 ， 而 且 积 分 亮度 的 质心 应 该 等 于 电子 面 质量 分 布 的 重心 在 双 极 
面 上 的 投影 ， 然 后 根据 TP 重 构 得 到 每 个 相交 极 面 上 的 CME 质心 分 布 ， 并 研究 CME 质心 


I 的 传播 方向 。 

= #2 Mierla 等 人 对 2007 年 5 月 15 日 事件 使 用 三 种 重 构 方法 获得 的 结果 

期 ” ”仪器 方法 平均 经 度 /(*) “平均 纬度 /(?) 
2007-05-15 COR1 LCT(from A)-TP 一 68 土 3 15 
2007-05-15 COR1 LCT(from B)-TP 一 73 士 3 15 
2007-05-15 CORI CM-TP 一 65 士 3 7 
2007-05-15 CORI PR-A 一 72 8 
2007-05-15 CORI PR-B 一 74 10 
2007-05-15 COR2 LCT(from A)-TP 一 61 士 3 11 
2007-05-15 COR2 LCT(from B)-TP 一 65 士 3 13 
2007-05-15 COR2 CM-TP 一 66 士 3 14 
2007-05-15 COR2 PR-A 一 53 13 
2007-05-15 COR2 PR-B 一 64 16 


Mierla 等 人 用 上 述 三 种 方法 获得 的 CME 纬度 和 经 度 的 结果 基本 一 致 ， 偏 差 不 超 过 
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10° (如 表 2 所 示 )。 但 不 同方 法 重 构 得 的 CME 传播 方向 (—51°, 79) 与 源 区 位 置 差异 较 大 ， 
这 可 能 是 由 于 CME 在 传播 过 程 中 发 生 了 径 向 偏转 所 导致 。 
4.2 2007 年 5 月 20 日 事件 
第 二 个 事件 是 发 生 于 2007 年 5 月 20 日 的 一 个 与 暗 条 爆发 相关 的 部 分 晕 状 CME, 
STEREO 两 个 航天 器 的 分 离 角 度 约 为 9%9。SOHO-LASCO 观测 得 到 的 CME 前 沿 的 投影 
度 为 275 km/s， 这 个 CME 的 传播 方向 指向 地 球 。Srivastava 45 AU 应 用 TP 方法 对 此 次 
事件 在 COR1 和 COR2 视 场 中 的 前 沿 进行 重建 ， 其 中 CME JD BJE2'S28 ， 前 沿 的 
速度 是 510 km/s， 而 重 构 得 到 的 天 空 平面 的 投影 速度 约 是 272 km/s. Mierla 等 人 ”利用 
3D-HT 技术 应 用 于 相同 的 特征 ， 重 建 结果 得 到 CME 源 区 位 置 为 B2°*S27*， 前 沿 的 三 维 速度 
是 548 km/s， 在 天 空 平面 的 投影 速度 约 250 km/s. 
TP 方法 与 3D-HT 方法 相似 ， 且 重 构 位 置 也 相同 ， 得 到 的 前 沿 的 三 维 速度 相近 。 
Srivasta-va ^& A" 在 重建 的 基础 上 进一步 计算 了 CME 的 三 维 传播 速度 以 及 到 达 地 球 的 时 
间 ， 这 与 地 球 附近 测 得 的 CME 实际 到 达 时 间 非 常 接近 ， 其 中 的 差异 可 能 是 由 测量 误差 以 及 
= 重 构 方法 存在 的 不 完备 性 所 引起 的 。 这 也 表示 如 果 简 单 使 用 投影 速度 估计 传播 时 间 的 话 会 
2 产生 很 大 的 误差 。 因 此 ， 对 CME 的 三 维 重 构 对 灾害 性 空间 天 气 预 报 具有 重要 意义 
4.3 2007 年 11 月 16 日 事件 
第 三 个 事件 是 2007 年 11 月 16 日 观测 到 的 一 个 西南 方向 的 CME， 当 时 STEREO 两 个 
航天 器 的 分 离 角度 为 40?。 这 个 事件 的 CME 具有 典型 的 三 分 量 结构 ， 即 包括 亮 前 沿 、 上 暗 腔 
和 亮 核 ，SOHO-LASCO 得 到 的 CME 前 沿 的 投影 速度 为 326 km/s. 
Howard 和 Tappin ”利用 TP 方法 重 构 了 该 事件 中 CME 前 沿 的 中 间 部 分 ， 得 到 的 结 
果 给 出 投影 速度 为 274.11 km/s， 三 维 速度 为 322.71 km/s, CME 源 区 位 置 是 S14°W73°. 
Temmer 等 人 四 分 别 利 用 STEREO-A+SOHO 和 STEREO-B+SOHO 的 日 冕 仪 数据 对 此 
事件 的 前 沿 进行 三 角 测 量 重 构 。 根 据 两 组 数据 反 演 得 到 CME 源 区 位 置 的 平均 值 为 
W120°S10?*， 三 维 速度 的 平均 值 为 403 km/s. Liu 等 人 ”根据 他 们 的 几何 三 角 测量 方法 
究 了 从 太阳 到 1 AU 的 CME 运动 学 (传播 方向 和 径 向 距离 )， 将 CME 成 像 重 构 与 MC 就 
ji. 地 重 构 结 果 进 行 对 比 ， 并 计算 出 CME 的 源 区 位 置 为 87"W123"，CME 速度 为 388 km/s. 
= Thernisien 4& A? 3x Jl GOS 方法 对 此 事件 的 磁 通 量 绳 进行 重 构 得 到 CME 的 源 区 位 置 是 
"T Sl4eW123。， 三 维 速度 为 345 km/s， 投 影 速度 为 289 km/s. 
上 述 结果 表明 几 种 重 构 方 法 得 到 的 速度 都 在 350 ~ 400 km/s 范围 内 ， 三 维 速度 都 略 高 
于 投影 速度 ， 但 各 自得 到 的 速度 都 不 完全 相同 。 这 种 差异 主要 是 因为 三 种 方法 重 构 针对 的 
CME 位 置 有 所 不 同 。 
4.4 2008 年 3 月 25 日 事件 
第 四 个 事件 是 2008 4E 3 月 25 日 的 一 次 与 暗 条 爆发 有 关 的 CME， 此 时 ，STEREO 两 个 
航天 器 的 分 离 角 度 是 47?，LASCO 给 出 的 CME 前 沿 的 投影 速度 为 1103 km/s， 源 区 位 置 为 
S10°E86°。 
Liewer 等 人 "利用 TP 方法 对 这 次 事件 进行 重 构 ， 得 到 前 沿 的 三 维 速度 是 1087 km/s, 
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前 沿 位置 S9°E86°. Mierla 等 人 四 


于 COR1 数据 ， 利 用 LCT-TP 方法 和 CM-TP 获得 的 源 区 
而 对 于 COR2 数据 两 种 方法 
其 它 方法 在 经 度 方面 差别 高 达 40"， 这 可 
射 的 偏振 度 过 低 ， 基 于 目标 4 
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将 LCT-TP, M-TP 和 PR 三 种 方法 应 用 到 这 个 事件 。 


的 结果 分 别 是 $S2 °E88°, S7°E92°. PR 方法 获得 的 源 
为 CME 中 心 
ij 振 状态 的 PR 重 构 方法 的 可 靠 
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对 
位 置 分 别 是 S15°E89°, S14°E88°, 

区 位 置 与 
高 密度 等 离子 体 ( 暗 条 物质 ) 辐 
二 程度 自然 会 受到 影响 。Mierla 等 


人 的 主要 目的 是 研究 CME 的 传播 角度 ， 所 以 并 未 给 出 
到 CME 的 速度 是 1130 km/s， 源 


等 人 四 利用 GCS 方法 对 CME 的 磁 通 量 绳 进行 重 构 ， 得 
区 位 置 S12°E84°。 

不 同方 法 对 这 个 事件 的 重 构 结果 给 出 
相差 不 超过 10。， 在 误差 范围 内 。Liewer 等 人 " 和 Thern 


Ww 


重 构 结 构 的 三 维 速度 。Thernisien 


的 源 区 位 置 基本 一 致 ， 重 构 结 果 与 实际 观测 结果 


isien 等 人 使 用 的 方法 虽然 差别 


明显 ， 但 给 出 的 CME 前 沿 位 置 的 速度 差别 不 大 ， 大 约 都 是 1100 km/s。 这 几 种 重 构 方法 得 


到 的 CME 传播 方向 偏差 都 在 10° 以 内 ， 而 Mierla 等 人 ™ 
到 的 传播 方向 相差 12° 左 右 ， 这 可 能 是 因 
4.5 2010 年 8 月 7 日 事件 
第 五 个 事件 是 2010 年 8 
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的 一 个 晕 状 CME, Ji 


利用 COR1 数据 和 COR2 数据 得 


为 CME 在 源 区 附近 发 生 了 偏转 。 


If, STEREO 两 个 航天 器 的 分 离 


角度 大 约 是 150*，LASCO 给 出 的 CME 前 沿 的 投影 速度 为 871 km/s. 


Feng 等 人 "将 几 种 不 同 
GCS 方法 以 及 局 部 相关 跟踪 加 三 角 测 
用 GCS 方法 获得 的 结果 较 准 确 
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(LCT-TR) 方法 。 


的 重建 方法 应 用 到 这 次 事件 上 


o 这 些 方法 分 别 是 MF 方法 、 
根据 已 有 的 结果 与 积累 的 经 验 ， 


， 因 此 本 文 以 GCS 方法 为 标准 来 验证 MF 方法 的 可 靠 性 ， 同 


时 计算 了 利用 不 同 重 构 方法 获得 的 三 维 CME 重心 及 其 经 红 


l 度 延展 范围 (如 表 3 所 示 )。 


#3 X CME 进行 重 构 得 到 的 相应 参数 /(°) 
方法 经 度 ”纬度 “最 小 经 度 ”最 大 经 度 “最 小 纬度 ”最 大 纬度 
MF 一 13.3 ”一 7.7 一 48 一 24.7 一 42.5 25.35 
GCS 13.4 7.3 35.6 —11 —42 27.6 
LCT-TR 27.8 2.8 —29.7 89.4 一 39.2 39.2 


由 表 3 可 以 看 出 ，MEF 方法 得 到 的 CME 重心 经 纬度 与 GCS 方法 得 到 的 相近 ，CME 延 


展 的 纬度 范围 相差 不 大 ， 但 是 经 度 延 展 的 范围 差别 明显 。 
偏差 都 很 大 ， 这 是 


AE 


这 种 相关 性 在 STEREO 双星 分 离 角度 过 大 时 的 可 靠 性 较 人 


因为 LCT 方法 的 基础 是 两 幅 STEREO 图 


LCT-TR 方法 得 到 的 所 有 结果 的 
像 上 等 离子 体 波动 的 相关 性 ， 
17 日 事件 发 生 时 ， 


氏 。2010 年 8 


STEREO-A 和 STEREO-B 之 间 的 分 离 角 度 约 150"， 导 致 LCT-TR 方法 重 构 的 结果 出 现 较 


大 偏差 。 
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CME 是 行星 际 环境 和 空间 天 气 扰动 的 主要 源头 ， 对 灾害 性 空间 天 气 的 预报 需要 实时 预 
测 CME 或 ICME 的 传播 方向 和 到 达 地 球 的 时 间 。 虽 然 ICME 到 达 地 球 的 时 间 是 灾害 性 空间 
天 气 预 报 的 重要 内 容 ， 但 是 更 重要 的 是 ICME 到 达 地 球 时 的 速度 和 南 向 磁场 分 量 (Bs); n 
外 ， 有 关 ICME 的 三 维 结构 和 三 维 速度 的 信息 也 具有 重要 的 参考 价值 。 本 文 详细 介绍 了 目 
前 常用 的 几 种 CME 成 像 重 构 方 法 和 两 种 CME 驱动 激 波 的 重 构 方 法 。 这 些 方法 各 自 以 不 同 
的 角度 强调 和 重建 CME 的 某 些 重要 特征 ， 但 是 对 其 他 方向 的 描述 则 与 实际 情况 相差 较 大 。 
这 说 明 以 现 有 的 观测 手段 获得 的 信息 还 不 能 反映 CME 的 全 貌 ， 也 难以 从 中 提取 出 CME 的 
完整 信息 。 新 的 探测 方法 、 观 测 手段 值得 进一步 探测 和 研发 ! 
目前 使 用 的 CME 成 像 重 构 方 法 分 为 两 个 大 类 : 基于 中 小 视 场 日 冕 仪 观测 数据 的 重 构 
和 基于 大 视 场 或 超大 视 场 日 冕 仪 (或 日 球 成 像 仪 ) 观测 数据 的 重 构 。 而 根据 不 同 的 数学 方 
法 和 所 针对 的 物理 特征 ， 各 个 大 类 中 又 包含 了 不 同 的 小 类 。 本 文 介绍 的 基于 日 冕 仪 数据 重 
构 CME 的 方法 包括 圆锥 模型 法 、 偏 振 比方 法 、GCS 模型 法 、3DCORE 模型 法 、 三 角 测 量 
ik. TP 方法 和 MF 方法 ; 除了 偏振 比方 法 之 外 ， 其 他 方法 都 可 以 得 到 CME 的 三 维 几何 结 
构 。 其 中 圆锥 模型 是 一 种 早期 描述 CME 几何 结构 的 模型 ， 它 没有 假设 CME 的 前 沿 特征 ， 
省 且 更 适合 于 描述 晕 状 CME; 但 是 圆锥 模型 中 没有 磁场 结构 ， 因 此 ， 该 方法 难以 可 靠 地 预 
ill CME 在 地 球 附 近 或 其 他 空间 位 置 处 的 磁场 ;， 不过， 这 个 模型 在 预测 CME 驱动 激 波 是 否 
到 达 地 球 以 及 到 达 时 间 方 面 获得 了 成 功 。GCS 对 于 单个 CME 的 三 维 速度 和 传播 方向 的 重 
建 效果 很 好 ; 但 是 在 针对 CME 连续 爆发 的 情况 时 ， 拟 合 精度 很 低 。 利 用 3DCORE 模型 进 
行 重建 时 ， 在 初始 参数 的 设置 中 会 经 常 使 用 GCS 的 拟 合 结果 ， 重 建 结果 主 要 是 1 AU 处 的 
磁场 信息 ， 这 对 地 磁 暴 的 预测 尤为 重要 。 偏 振 比 方法 是 将 一 个 视点 的 偏振 观测 转化 为 与 天 
空 平 面 的 3D 距离 ， 只 能 得 到 沿 每 条 视线 的 CME 等 离子 体 密度 的 距离 加 权 平 均 ， 不 能 提供 
CME 的 深度 信息 ， 也 不 能 提供 CME 的 速度 和 位 置信 息 。 三 角 测 量 和 TP 方法 不 需要 假设 
CME 的 几何 结构 ， 只 要 CME 的 特征 在 图 像 中 可 以 清晰 地 分 辨 ， 就 可 以 得 到 比较 好 的 结果 ， 
但 是 这 两 种 方法 都 不 能 重 构 CME 中 的 电子 密度 。 而 且 TP 方法 的 一 个 假设 是 仿 射 几何 的 基 
本 框架 成 立 ， 所 以 会 有 一 个 放大 因子 ， 得 到 的 结果 与 实际 情况 会 有 一 定 偏 差 。 由 TP 方法 
发 展 的 HT 技术 有 一 个 独特 的 优点 ， 可 以 即时 显示 运动 过 程 及 其 速度 并 且 灵 敏 度 高 ， 可 以 
提取 最 模糊 的 特征 并 检测 它们 的 运动 。MEF 方法 在 不 假设 CME 形状 的 情况 下 ， 能 很 好 地 重 
建 CME 表面 ， 使 用 起 来 更 加 灵活 ;在 获得 CME 三 维 形状 之 后 可 以 分 析 其 几何 中 心 与 主 旬 
尺度 等 ， 并 且 MF 方法 是 基于 三 个 视角 进行 重建 的 ， 比 基于 两 个 视角 的 重建 得 到 的 结果 更 
准确 ，MF 方法 的 缺点 是 所 得 的 结果 没有 包括 CME 的 内 部 结构 和 等 离子 体 分 布 。 
基于 日 球 成 像 仪 数据 的 重 构 方法 包括 几何 三 角 测 量 、SSE 方法 、 椭 圆 演 化 模型 法 和 相 
关 辅 助 重 构 法 。SSE 方法 和 椭圆 演化 模型 都 要 假设 CME 的 几何 结构 。SSE 方法 在 CME 4E 
展 角 过 大 或 者 过 小 的 时 候 都 会 导致 结果 不 可 信 ; 而 延展 角 不 能 过 大 将 限制 SSE 模型 用 于 空 
间 天 气 预报 ;， 男 外，SSE 模型 假设 CME 的 前 沿 是 圆 形 截 面 。 顶 圆 演化 模型 假设 截面 为 椭 医 
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这 种 修正 使 得 模型 更 灵活 ， 更 适合 多 点 观测 的 自治 建 模 。 几 何 三 角 测量 相对 于 前 面 这 两 

种 方法 的 独特 性 和 优点 是 : 没有 自由 参数 ， 所 做 假设 少 ， 并 且 可 以 对 CME 从 太阳 连续 追踪 

至 1 AU。 相 关 辅助 重 构 法 则 针对 小 尺度 太阳 风 瞬 变 所 提出 ， 对 于 大 尺度 CME 的 重 构 则 是 

将 其 分 解 为 无 数 个 小 CME。 这 个 方法 的 优点 是 不 需要 假设 CME 的 几何 结构 并 可 以 自动 化 

进行 ， 操 作 简单 ， 还 可 以 得 到 三 维 空间 中 OME 的 径 向 速度 分 布 ， 缺 点 是 只 能 应 用 于 日 地 连 

线 附近 的 CME， 且 共 面 效应 会 导致 共 面 附近 CME 的 纵向 定位 精度 非常 低 。 

CME 驱动 激 波 的 三 维 重 构 也 非常 重要 ， 因 为 CME 与 其 激 波 紧 密 相关 ， 而 激 波 重 构 可 
以 给 出 激 波 的 三 维 结构 和 运动 学 。 激 波 椭 球 模型 大 加 GCS 模型 再 应 用 三 个 不 同 视角 的 观测 
数据 ， 可 以 很 好 地 再 现 CME 及 其 驱动 激 波 的 三 维 结构 。 激 波 球体 模型 可 应 用 到 极端 CME 
产生 的 激 波 ， 这 种 CME 速度 非常 快 ， 驱 动 的 激 波 角 宽度 非常 宽 。 目 前 ， 将 激 波 重 构 结 果 与 
OME 特性 、 高 能 粒子 特性 、 太 阳 风 局 地 观测 作对 比 ， 对 空间 天 气 的 变化 和 高 能 粒子 加 速 的 
研究 都 具有 重要 意义 ， 也 是 空间 科学 和 太阳 物理 研究 中 备 受 关注 的 热门 课题 。 

如 前 所 述 ， 无 论 哪 种 拟 合 都 含有 多 个 自由 参数 ， 这 些 自由 参数 之 间 是 简 并 的 (特别 是 非 

J 线性 拟 合 )， 所 以 解 不 唯一 ， 也 不 稳定 ， 自 由 参数 越 多 ， 解 越 不 确定 。 因 此 ， 引 入 更 多 自由 

参数 ， 虽 然 带 来 一 定 的 自由 度 ， 但 方法 的 缺点 也 更 明显 。 这 也 是 没有 自由 参数 的 几何 三 角 测 

量 方法 的 优越 性 所 在 。 这 些 重 构 方法 在 对 CME 的 运动 学 、 动 力学 的 研究 方面 发 挥 了 重要 的 

作用 ， 但 它们 都 存在 一 些 缺点 。 首 先 ， 这 些 重 构 方法 都 做 了 很 多 假设 并 加 上 了 相应 的 约束 条 

fb. 其 次 ， 根 据 三 维 图 的 结构 原理 ， 由 二 维 图 像 构造 一 个 目标 的 三 维 结构 至 少 需要 三 幅 该 目 

标的 二 维 图 像 ， 而 且 获 得 这 三 幅 二 维 图 像 的 观测 点 与 目标 的 连 线 不 能 位 于 同一 平面 上 ， 并 且 

这 三 条 连 线 最 好 是 相互 垂直 。 但 是 目前 我 们 对 太阳 及 其 活动 的 成 像 观 测 基本 上 都 是 在 黄道 

面 上 进行 的 ， 所 以 依据 这 样 的 观测 数据 对 我 们 所 研究 的 目标 进行 三 维 重 构 时 ， 显 然 还 缺少 一 

个 维度 的 信息 。 

CME 内 部 磁场 结构 、 密 度 极 值 、 磁 场 B. 分 量 即 南 向 磁场 分 量 都 是 研究 CME 物理 本 

质 的 重要 参数 ， 其 中 南 向 磁场 分 量 对 预测 地 磁 暴 是 一 个 非常 关键 的 量 。 而 CME 成 像 重 构 一 

般 无 法 提供 OME 内 部 的 磁场 信息 ， 所 以 须 先 利用 一 些 基于 局 地 观测 数据 的 磁 云 重 构 方 法 来 

计算 磁场 相关 信息 。 比 如 最 早 的 柱 对 称 无 力 场 磁 通 量 绳 模型 一 ， 以 及 近期 的 椭圆 -圆柱 

解析 磁 通 量 绳 模型 "都 可 以 给 出 磁 云 中 磁场 各 个 分 量 以 及 其 他 相关 信息 。 还 有 一 种 不 需要 
假设 磁 云 横 截 面 的 Grad-Shafranov (GS) 方法 ， 起 初 是 用 于 研究 地 球 磁 层 的 磁场 结构 

后 来 Hu 等 人 一 以 及 Liu 等 全 "发现 该 方法 可 以 用 于 磁 云 中 的 磁场 分 析 。 而 GS 3 

构 的 结果 则 可 以 通过 具有 分 离 角度 的 多 个 航天 器 观测 数据 进行 验证 ， 这 要 求 几 颗 卫 星 分 离 

距离 较 远 (但 也 不 能 太 远 ， 否 则 MC 会 被 太阳 风 扭 曲 ) 才能 给 出 观测 的 一 致 性 ， 因 此 Liu 等 

AE STEREO 和 Wind/ACE 几 颗 卫星 的 观测 结果 进行 验证 。 以 前 验证 只 用 Wind 和 

ACE， 由 于 这 两 个 卫星 距离 太 近 ， 本 质 上 观测 到 的 是 同一 结构 ， 无 法 进行 有 效 验 证 。 还 有 

一 些 与 磁 云 重 构 相 关 的 早期 代表 性 工作 ， 可 以 参见 文献 ID OCHO). SR, CME 成 像 重 

构 与 磁 云 (或 ICME) 局 地 重 构 之 间 的 对 比 研究 将 受到 越 来 越 多 的 关注 ， 因 为 这 可 以 给 出 磁 

通 量 绳 结构 和 磁场 定向 的 重要 信息 ， 如 磁 通 量 绳 轴 的 定向 、 在 运动 过 程 中 是 否 发 生 了 旋转 

等 。 目 前 ， 一 些 学 者 已 经 开展 这 方面 的 工作 ， 如 Liu 等 人 利用 GCS 成 像 重 构 获得 CME 
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的 整体 结构 ， 再 利用 GS 局 地 重 构 方法 获得 1 AU 附近 局 部 磁 通 量 绳 的 相关 信息 ;， Chen 等 
A 也 有 类 似 的 工作 等 。 灾 害 性 空间 天 气 的 预报 不 仅 需 要 CME 到 达 地 球 的 时 间 预 报 ， 更 
需要 南 向 磁场 分 量 (B) 和 速度 的 预报 。 这 就 体现 了 CME 成 像 重 构 结 果 和 磁 云 (或 ICME) 
局 地 重 构 结果 之 间 对 比 的 重要 性 。 

获取 这 些 参 数 的 正确 数值 有 赖 于 对 CME 三 维 结构 的 完备 和 正确 的 重 构 。 完 备 的 意思 
是 指 在 重 构 的 过 程 中 不 需要 加 入 任何 假设 ， 也 不 需要 依赖 任何 模型 ， 只 依赖 于 来 自 三 个 相 
互 垂 直 的 视角 的 观测 资料 。 当 然 ， 目 前 的 观测 资料 还 不 能 满足 这 个 要 求 ， 因 此 ， 所 使 用 的 
重 构 方法 都 包含 相应 的 假设 或 者 依赖 于 某 个 模型 。 克 服 这 些 困难 的 唯一 出 路 在 于 开展 在 黄 
道 面 以 外 的 观测 ， 而 2020 年 欧洲 发 射 的 太阳 轨道 器 (solar orbiter) T 则 实现 了 黄道 面 外 的 
观测 。 太 阳 轨 道 器 携带 日 园 仪 和 日 球 仪 与 太阳 的 距离 约 0.28 AU， 虽 然 没 有 比 帕克 太阳 探 
针 (Parker Solar Probe) ^ 更 接近 太阳 ， 但 它 相 对 于 黄道 面 的 倾角 更 大 ， 能 够 从 更 高 纬度 
对 太阳 进行 观测 。 所 以 利用 太阳 轨道 器 的 观测 数据 重 构 三 维 CME， 将 会 为 我 们 揭示 更 多 的 
CME 三 维 信息 ， 帮 助 我 们 更 加 深入 、 全 面 地 研究 CME 的 内 部 磁场 结构 和 等 离子 体 分 布 。 
如 果 还 能 够 对 太阳 进行 抵 近 观 GR) W, RIII CME/ICME 的 重 构 过 程 就 会 更 可 靠 ， 据 
此 进行 的 相应 预报 也 将 会 更 加 准确 。 
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Advances in 3D Reconstruction of Coronal Mass Ejections 


ZHAO Xing-mei!?, FENG Li, SONG Hong-qiang +, LIN Jun!? 


(1. Yunnan Astronomical Observatory, Chinese Academy of Sciences, Kunming 650011, China; 2. 
University of Chinese Academy of Sciences, Beijing 100049, China; 3. Purple Mountain Observatory, 
Nanjing 210023, China; 4. Shan Dong University, Weihai 264209, China) 


Abstract: Coronal mass ejection (CME) is the large scale magnetized plasmoid ejected 
from the Sun, which brings huge amount of magnetic flux and plasma into interplanetary 
space. An earthward CME will interact with the magnetosphere of the Earth and invokes the 
substorm and the other phenomena of the space weather as it approaches to the Earth. The 
2-dimansional data provided by the current observational techniques cannot describe the true 
magnetic structure and the plasma distribution of CMEs comprehensively. We need to look 
into the 3-dimensional structure and the associated three components of CME speeds in order 
to predict the time when an ICME reaches the Earth and the potential consequent impact on 
the Earth and the nearby environment. In this paper, 3D reconstruction methods of CME 
based on existing imaging observations are introduced, including two kinds of reconstruction 
methods based on coronagraph data and heliosphere imager data and CME-driven shock 
wave 3D reconstruction methods with high correlation with CME imaging reconstruction. 
Each method shows apparent advantages in dealing with specific events, but its weakness 
and necessary constrains to its applications exist as well. Results obtained via various 
methods are compared in this work, and we found that CME velocities and moving directions 
deduced from these methods are fairly close to one another, which shows high reliability of 
these methods. Finally, the hot topics related to the 3-dimanesional reconstruction of CME 


(ICME) and the relevant development in reconstructing methods are also discussed. 


Key words: corona; coronal mass ejection; 3D reconstruction 


